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SUMMARY 


A  method  is  developed  for  predicting  the  nonperiodic  air  loads  caused  by 
control  inputs  applied  to  a  rotary  wing  in  forward  flight.  Use  is  made 
of  a  numerical  description  of  the  geometry  and  circulation  strength  of  the 
vorticity  in  the  wake  to  compute  the  time -varying,  nonuniform -flow 
field  in  the  plane  of  the  rotor  disc.  The  method  of  approach  is  similar 
to  that  developed  previously  for  the  steady-state  flight  condition,  the 
major  difference  being  that  the  solution  procedure  lends  itself  to  the 
treatment  of  transient  phenomena,  such  as  the  nonperiodic  loadings 
caused  by  time-varying  collective  pitch. 

Several  approximations  are  made  to  simplify  the  analysis.  The  major 
assumptions  are:  (1)  the  rotor  blades  are  structurally  rigid  and  have 
only  a  flapping  degree  of  freedom;  (2)  the  hub  of  the  rotor  continues  to 
translate  in  level,  constant-speed  flight  during  the  short  time  interval 
of  interest;  (3)  the  geometry  of  the  wake  can  be  specified  a  priori;  and 
(4)  no  account  need  be  taken  of  the  shed  vorticity  in  the  wake. 

Computed  flapping  and  air  load  distributions  are  compared  with  tran¬ 
sient  data  obtained  in  wind-tunnel  tests  on  a  full-scale  H-34  rotor.  In 
general,  very  encouraging  agreement  is  found  leading  to  the  conclusion 
that  calculation  and  prediction  of  nonperiodic  loadings  on  rotary  wings  is 
a  feasible  task. 


FOREWORD 


This  investigation  of  air  loadings  on  rotary  wings  was  conducted  at  the 
Cornell  Aeronautical  Laboratory,  Inc.  (CAL),  under  U.  S.  Army 
Contract  DA  44-177-AMC-77(T)  during  the  period  September  1963 
through  May  1965.  The  program  was  sponsored  by  the  U.  S.  Army 
Aviation  Materiel  Laboratories  (USAAVLABS),  Fort  Eustis,  Virginia, 
as  Task  1P125901A142,  and  was  administered  by  Mr.  J.  E.  Yeates. 

Mr.  Leonard  Segel  is  the  project  engineer  and  author  of  this  report. 

He  was  aided  in  a  very  large  measure  by  Mr.  Edwin  F.  Chmielewski  of 
the  Computer  Services  Department  (CAL)  who  performed  the  mon¬ 
umental  task  of  generating  a  working  computing-machine  program. 
Acknowledgement  is  also  made  of  the  combined  efforts  of  personnel 
from  the  Ames  Research  Center,  National  Aeronautics  and  Space 
Administration  and  Sikorsky  Aircraft  Division,  United  Aircraft 
Corporation,  in  providing  the  experimental  data  used  to  validate  the 
analysis. 


v 


CONTENTS 


Page 

SUMMARY  iii 

FOREWORD  v 

LIST  OF  ILLUSTRATIONS  ix 

LIST  OF  SYMBOLS  xi 

INTRODUCTION  1 

CONCLUSIONS  3 

RECOMMENDATIONS  5 

PRIOR  STUDIES  OF  TRANSIENT  AIR  LOADS  ON  ROTARY 

WINGS  7 

REPRESENTATION  OF  THE  ROTOR  WAKE  10 

GENERAL  REMARKS  10 

DEVELOPMENT  OF  A  SIMPLIFIED  WAKE  MODEL  13 

THE  COMPUTATIONAL  MODEL  16 

COMPARISON  OF  THEORY  AND  EXPERIMENT  18 

STEADY -STATE  LOADINGS  18 

RESPONSE  TO  CHANGE  IN  COLLECTIVE  PITCH  19 

ADDITIONAL  NUMERICAL  INVESTIGATIONS  23 

REFERENCES  54 

APPENDIX-  THE  MATHEMATICAL  MODEL  AND  57 

COMPUTATIONAL  PROCEDURE  USED  TO 
CALCULATE  TRANSIENT  BLADE  LOADS 

Indexing  Nomenclature  57 

vii 


Input  Data  57 

Wake-Geometry  Computations  58 

Computation  of  <r~  Matrix  51 

Induced -Velocity  Computation  62 

Rotor  Kinematics  63 

Solution  Procedure  for  Unknown  Circulation  Strengths  65 
Blade -Flapping  Dynamics  68 

Digital-Computer  Program  69 

DISTRIBUTION  71 


viii 


ILLUSTRATIONS 


Figure  Page 


1  Pictorial  Example  of  the  Initial  Portion  of  the  Wake  28 

of  a  Two-Bladed  Rotor  Divided  Into  Four  Radial 
Segments. 

2  Wake  Models  Examined  in  a  Wake -Simplification  28 

Study. 

3  Influence  of  Wake  Representation  on  Computed  29 

Aerodynamic  Loading  vs.  Azimuth  Angle. 

4  Comparison  of  Calculated  and  Measured  30 

Loadings  in  Hovering  Flight. 

5  Schematic  Diagram  of  TBL  Model.  31 

6  Axis  System  Used  in  TBL  Analysis.  32 

7  Steady-State  Aerodynamic  Loadings  vs.  33 

Azimuth  Angle;  H-34  Rotor,  >f=  0. 18. 

8  Comparison  of  Induced -Velocity  Distributions  34 

Yielded  by  SSBL  and  TBL  Calculations* 

9  Steady-State  Aerodynamic  Loadings  vs.  35 

Azimuth  Angle, 

10  Measured  and  Calculated  Transient  Flapping  and  36 

Air  Load  Responses;  H-34  Rotor,  Wind-Tunnel  Run 

No.  l,/4s  0.  202. 

11  Measured  and  Calculated  Transient  Flapping  and  37 

Air  Lead  Responses;  H-34  Rotor,  Wind-Tunnel  Run 

No.  2./*  =  0.  286. 

12  Measured  and  Calculated  Transient  Flapping  and  38 

Air  Load  Responses;  H-34  Rotor,  Wind-Tunnel  Run 

No.  3,  ^4  —  0.287. 


ix 


13  Total  Lift  on  Blade  vs.  Azimuth;  H-34  39 

Rotor(/44  =  0.18. 

14  Induced  Velocity  vs.  Span;  Computed  for  Two  40 

Different  Distributions  of  Vorticity  in  the  Wake; 

H-34  Rotor, =  0. 18. 

15  Transient  Response  to  a  Rapid  Increase  in  Collective  41 

Pitch;  Calculated  for  the  H-34  Rotor  in  Hovering 
Flight. 

16  Transient  Response  to  a  Rapid  Decrease  in  Collective  42 


Pitch;  Calculated  for  the  H-34  Rotor  in  Hovering 
Flight. 

17  Indexing  Procedure  Used  to  Designate  Blade  and  43 

Wake  Geometry. 

18  Condensed  Flow  Charts:  Main  Program  -  "Blades".  44 

19  Condensed  Flow  Charts:  Subroutine  -  "Comequ".  48 

20  Condensed  Flow  Chart:  Subroutine  -  "Proced",  51 

21  Condensed  Flow  Chart:  Subroutine  -  "Comthe".  52 

22  Condensed  Flow  Chart:  Subroutine  -  "Comwb".  53 


x 


LIST  OF  SYMBOLS 


</J<? 

LCA/ 

HBLM 

A/A 

A/S 

A/// 

R 

T 

Or 

Or 


Blade  segment  number  for  maximum  circulation 
Moment  of  inertia  of  blade  about  flap  hinge,  slug -ft  ^ 
Number  of  spanwise  blade  segments 

Moment  about  flapping  hinge  caused  by  aerodynamic  loading, 
ft /lbs 

Azimuth  index  limit  for  nonrolled  wake 

Weight  moment  of  blade,  ft  /lbs 

Number  of  azimuth  positions  per  rotor  revolution 
(an  integral  multiple  of  A/B ) 

Number  of  blades 

Number  of  revolutions  of  wake 


Rotor  radius,  ft 


Aerodynamic  loading  at  blade  load  point,  lbs/ ft 


Component  of  air  velocity  relative  to  the  blade  (at  the 
three-quarter  chord  position)  in  a  plane  containing  the 
blade -span  axis  and  the  rotor  axis  and  oriented  perpendicular 
to  the  blade  span,  ft/sec 

Component  of  air  velocity  relative  to  the  blade,  perpendicular 
to  the  plane  containing  the  blade-span  axis  and  the  rotor  axis, 
ft/ sec 


xi 


V  Forward  flight  velocity,  ft/ sec 

X,Y,Z  Shaft  axia  system  with  origin  located  at  rotor  hub 
(see  Figure  6) 

nondimensional 

b  Blade  semichord,  ft 

e  Radial  distance  from  rotor  axis  to  blade  root,  ft 

q  Gravitational  acceleration,  ft /sec  ^ 

h  "Effective  plunging"  velocity  of  the  three-quarter  chord 
point  of  a  specified  segment  on  a  rotor  blade,  ft/ sec 

i  An  index  designating  a  blade  on  a  rotary  wing 

U  An  index  designating  the  wake  trailed  by  blade 

j  An  index  designating  a  spanwise  position  on  a  rotor  blade 

J  An  index  designating  a  spanwise  position  in  the  wake 

4**V*l’a  Components,  along  the  axes,  respectively,  of  the 

position  vector  between  the  "forward"  end  of  a  wake  element 
and  a  load  point  on  the  blade,  ft 

Components,  along  the  axes,  respectively,  of  the  position 
vector  between  the  "aft"  end  of  a  wake  element  and  a  load 
point  on  the  blade,  ft 

m  An  index  designating  an  instant  of  time' 


Lift  curve 


,lope  (&)  • 


xii 


An  index  designating  an  azimuthal  position  in  the  wake 


Radial  distance  from  center  of  rotation  to  midpoint  of  blade 
segment,  or  to  other  specified  point,  ft 


Radial  distance  from  rotor  axis  to  end  point  of  blade 
segment,  ft 


Radial  distance  from  rotor  axis  to  spanwise  location  of 
rolled -up  root  and  tip  vortex,  ft 


Length  of  position  vector  between  a  load  point  and  the 
forward  and  aft  ends,  respectively,  of  a  wake  element,  ft 


Nondimensional  radial  distance  from  rotor  axis  to  specified 
point  on  a  blade 


*,«/.  y  components  of  velocity  induced  by  the  wake  at  a  specified 
point  on  a  blade,  ft/ sec 


Mean  induced  downwash  governing  the  vertical  transport  of 
vorticity  in  the  wake,  ft/ sec 


y  component  of  velocity  induced  by  all  elements  of  the  wake 
other  nhan  the  trailing  elements  adjacent  to  the  blades 


A  "wind"  axis  system  with  the  origin  located  at  the  hub  of  the 
rotor  (see  Figure  6) 

Circulation,  either  bound  or  trailing,  ft  ^/ sec 


Angular  velocity  of  rotor  shaft,  rad/ sec 


Angle  of  attack  at  specified  blade  segment,  rad 


Shaft  angle  of  attack  (positive  for  forward  tilt  of  the  rotor  shaft), 
rad 


p  Flap  angle  of  blade  relative  to  plane  normal  to  rotor  shaft,  rad 


€  Flap  hinge  offset,  ft 

0  Pitch  angle  of  a  given  blade  segment  (as  measured  in  a  shaft 

axis  system),  rad 

&0  Collective  pitch  setting  of  rotor  (reference  is  at  the  rotor 
axis)  rad 


Cosine  cyclic  pitch,  rad 


0f$  Sine  cyclic  pitch,  rad 


Jd/ff  Blade  twist  (linear),  rad/ft 


Coordinates  of  endpoints  of  trailing  vorticity  as  measured  in  the 
it ,tf.y  axis  system,  ft. 

0~  Influence  coefficient  yielding  induced  velocity  at  a  specified  point 
as  caused  by  an  element  of  vorticity  of  unit  strength,  1/ft 


Angle  between  air  velocity  vector  relative  to  the  blade  and 
the  rotor  shaft  plane,  rad 


f  Azimuth  angle  (  p  =  0  when  blade  is  oriented  along  positive 
ft  axis),  rad 


NOTE:  A  single  dot  over  a  variable  denotes  the  first  time  derivative; 

a  double  dot  denotes  the  second  time  derivative. 


xiv 


INTRODUCTION 


The  design  of  rotor  blades  to  tolerate  the  oscillatory  stresses  encoun¬ 
tered  by  a  helicopter  in  forward  flight  is  recognized  as  one  of  the  major 
problems  facing  the  designers  of  helicopter  rotor  systems.  This 
problem,  together  with  the  requirement  for  improving  the  structural 
integrity  of  rotary  wings,  has  necessitated  the  development  of  research 
programs  whose  objective  is  the  derivation  of  rational  and  valid 
procedure  for  predicting  the  time-varying  air  loadings  and  stressing 
of  rotor  blades  designed  for  helicopter  applications.  Since  the 
problem  is  worldwide,  research  has  been  and  is  being  performed  in 
a  number  of  countries.  The  VTOL  Dynamic  Loads  Symposium 
held  in  June  1963  under  the  joint  sponsorship  of  the  Cornell  Aero¬ 
nautical  Laboratory,  Inc.  (CAL)  and  the  U.  S.  Army  Aviation  Materiel 
Laboratories  (USAAVLABS)  gave  evidence  of  the  effort  being  mounted 
both  in  the  United  States  and  abroad. 

It  is  now  recognized  that  the  crux  of  the  helicopter  blade-load  prediction 
problem  is  the  requirement  to  model  the  airflow  through  the  rotor 
in  considerably  more  accurate  detail  then  has  been  sufficient  hereto¬ 
fore  for  studies  of  rotor  performance.  To  tms  end,  investigators  such 
as  Willmer  [  1  ],  Molyneux[2],  Miller  [  3,  4,  5],  Tararine[6], 

Shi-Tsun  [7,  8],  Davenport  [9],  DuWaldt  and  Piziali  [10,  11]  have 
developed  mathematical  models  of  the  rotor -wake  system  (in  a 
variety  of  formats)  for  the  purpose  of  predicting  the  air  loadings 
produced  in  steady,  forward  flight.  It  appears  that  these  investiga¬ 
tions  can,  in  part,  be  differentiated  from  one  another  by  the  extent 
to  which  resort  is  had  to  high-speed  computing  machinery  in  order  to 
minimize  the  necessity  for  simplifying  the  mathematical  treatment. 

The  investigation  reported  herein  is  an  extension  of  the  rotor-wake 
modeling  concepts  previously  described  in  Reference  12.  This  extension 
consists  primarily  in  the  recasting  of  the  original  analysis  that  was 
suitable  for  steady- state  calculations  to  one  that  yields  a  temporal 
and  spatial  history  of  the  pertinent  response  variables  as  a  function  of 
time-varying  changes  in  collective  pitch.  Accordingly,  reference  shall 
be  made  throughout  this  report  to  "steady-state  blade  loads"  (SSBL) 
calculations  and  to  "transient  blade  loads"  (TBL)  calculations. 

It  should  be  understood  that  the  SSBL  calculation  (as  described  in 
Reference  12)  consists  of  a  procedure  wherein  the  instantaneous 
velocities  induced  by  the  wake  of  the  rotor  are  determined  for  each 
blade  section  in  a  rotating  coordinate  system.  Further,  the  SSBL  model 
of  wake  vorticity  results  in  an  inflow  distribution  that  varies  with  time 
and  position  in  the  rotor  disc,  but  which  is  periodic  with  each  revolution 
of  the  rotor  thus  producing  an  integrated  mean  value  of  inflow  that  does 
not  change  with  time.  Thus,  we  use  the  term  "steady  state",  even 
though  the  inflow  seen  by  a  blade  section  is  varying  continually  during 
one  revolution  of  the  rotor. 
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On  the  other  hand,  a  TBL  calculation  must  consist  of  a  procedure  where¬ 
in  nonperiodic  induced  velocities  and  loadings  can  be  determined.  When 
collective  blade-pitch  control  is  applied,  such  as  to  change  the  mean 
thrust  being  produced  by  a  rotor,  the  momentum  that  is  imparted  to  the 
air  mass  by  the  rotor  must  be  changed.  This  is  to  say  that  the  air  mass 
must  be  accelerated  until  a  new  steady-state  thrust  corresponding  to  the 
changed  operating  conditions  is  achieved.  From  the  point  of  view  of 
momentum  theory,  a  time  interval  must  elapse  during  which  the  air  mass 
associated  with  a  rotor  is  accelerated  to  its  new  steady  state  inflow 
velocity.  From  the  microscopic  point  of  view,  wherein  inflow  is  viewed 
as  being  induced  by  wake  vorticity,  additional  time  rates  of  change  of 
inflow  are  created  by  the  existence  of  a  wake  of  vortex  filaments  and 
strengths  that  are  no  longer  periodic  with  respect  to  one  revolution  of  the 
rotor.  Thus,  there  is  an  equivalence  between  the  momentum  concept  of 
an  air  mass  being  accelerated  during  the  transient  interval  between  two 
steady  states  and  the  concept  in  which  a  change  in  vorticity  level  is 
diffused  through  the  wake  bringing  about  a  change  in  the  velocities 
induced  at  the  rotor  disc  by  the  wake. 

Other  investigators  [13,  14]  have  measured  the  response  of  a  hovering 
or  vertically-moving  rotor  to  a  change  in  collective  pitch.  They  found 
that  the  thrust  response  was  highly  sensitive  to  the  rate  of  change  of 
collective  pitch  for  rapid  control  inputs  and  were  able  to  explain  the 
measured  results  in  terms  of  the  requirement  to  accelerate  an  apparent 
additional  mass  of  air  associated  with  a  rotor  considered  to  be  an 
impervious  disc.  Although  the  concept  of  additional  mass  may  suffice  to 
relate  the  time-varying  total  thrust  on  a  blade  or  rotor  to  a  time  variation 
in  blade  pitch,  this  concept  is  not  suitable  for  computing  the  inflow 
distribution  required  to  determine  blade  loadings  as  a  function  of 
spanwise  location,  azimuth  position,  and  time.  Rather  it  is  necessary 
to  compute  the  inflow  caused  by  a  wake  possessing  elements  of  vorticity 
whose  strengths  vary  in  a  nonperiodic  fashion  as  a  result  of  a  time- 
varying  change  in  collective  pitch.  It  is  this  particular  task  that 
consituted  the  major  objective  of  this  investigation.  Successful 
accomplishment  of  this  task  is  viewed  as  a  first  and  major  step  towards 
computing  the  detailed  temporal  and  spatial  distribution  of  aerodynamic 
loads  encountered  by  a  helicopter  rotor  during  maneuvering  flight. 

Earlier  analyses  of  transient  air  loads  on  rotary  wings  are  reviewed  in 
the  report  prior  to  presenting  the  analytical  approach  developed  for 
purposes  of  modeling  TBL  phenomena.  The  presentation  of  the  method  of 
approach  used  in  this  investigation  includes  a  discussion  of  rotor-wake 
representation  and  a  brief  description  of  the  developed  computational 
procedure.  This  description  is  followed  by  a  comparison  of  theoretical 
calculations  and  experiment,  with  the  technical  discussion  being  concluded 
by  presenting  both  the  current  capabilities  and  shortcomings  of  the  devel¬ 
oped  analytical  procedure,  as  applied  to  specific  flight  regimes.  Details 
pertaining  to  the  developed  mathematical  model  and  computational 
procedure  are  given  in  an  appendix. 
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CONCLUSIONS 


This  investigation  has  shown  that  the  prediction  of  nonperiodic  air 
loadings  on  rotary  wings  is  a  feasible  task.  Although  the  agreement 
achieved  between  theory  and  experiment  indicates  that  the  major  aspects 
of  transient  blade  loading  phenomena  are  adequately  represented  in  the 
analysis,  future  refinements  of  the  existing  mathematical  model  will  be 
required  to  provide  a  prediction  capability  having  a  high  order  of 
accuracy.  Ultimately,  it  should  prove  possible  to  use  the  TBL  calculation 
procedure  to  investigate  almost  any  rotor  dynamics  problem  wherein  it 
is  necessary  to  include  an  adequate  description  of  the  wake-induced 
flow  field.  This  work  should  be  looked  on  as  a  first  step  towards  the 
goal  of  computing  the  air  loadings  on  the  rotor  of  a  maneuvering  helicopter. 

It  is  further  concluded  that: 

(1)  A  wake  model  consisting  only  of  trailing  vortex  elements, 
trailed  as  a  consequence  of  a  spanwise  variation  in  bound 
vorticity,  constitutes  a  first-order  representation  of  the 
wake  of  a  rotary  wing. 

(2)  The  radial  contraction  of  the  rolled-up  root  and  tip 
vortices  of  a  rotary  wing  must  be  simulated  in  order 

to  obtain  accurate  predictions  of  the  air  load  distribution 
produced  in  hovering  flight.  It  would  appear  that  this 
conclusion  *lso  holds  for  forward  flight  at  very  small 
advance  ratios. 

(3)  For  the  limiting  case  of  steady  forward  flight,  the  TBL 
computational  procedure  provides  results  which  are 
equivalent  to  those  yielded  by  the  previously  developed 
SSBL  calculations.  Identical  steady  periodic  air  loadings 
are  predicted  by  both  procedures  provided  the  tip-path- 
plane  orientation  predicted  by  the  TBL  procedure  is  in 
agreement  with  the  orientation  specified  in  the  SSBL 
calculation.  As  a  result  of  the  integration  performed 

in  the  TBL  calculation  to  obtain  the  total  lift  on  the 
rotor,  it  was  observed  that  nonuniform -inflow  calculations 
yield  an  oscillatory  four -per -revolution  component  of 
total  shear  at  the  hub  having  an  amplitude  of  0.  1  g 
(approximately)  in  all  cases  studied  for  the  H-34, 
whereas  uniform -inflow  calculations  yield  very  little 
oscillatory  component  in  the  total  shear  transmitted  to 
the  rotor  hub. 

(4)  The  aerodynamic  model,  in  its  present  form,  together 
with  a  model  of  the  kinematics  and  dynamics  of  a  flapping 
degree  of  freedom,  yields  flapping  and  air  loading  responses 
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caused  by  a  change  in  collective  pitch  that  are  in  good 
agreement  with  wind  tunnel  measurements  made  on  a 
full-scale  H-34  rotor. 

(5)  Preliminary  and  tentative  support  has  hr  ^n  obtained  for  a 
hypothesis  stating  that  a  wake  assumed  t  have  a  constant 
strength  corresponding  to  the  mean  total  load  on  a  rotor 
would  produce  a  reasonable  approximation  of  the  velocity 
field  that  is  induced  in  the  plane  of  the  rotor  disc.  A 
considerable  number  of  computational  studies  is  required 
to  place  this  hypothesis  on  a  firmer  foundation. 

(6)  Although  the  wake  spacing  (i.  e.  ,  wake  pitch)  assumed  to 
exist  in  the  steady  state  preceding  and  following  a  transient 
interval  has  a  significant  influence  on  the  steady  total  load 
and  load  distribution  achieved  in  hover,  the  timewise 
variation  of  loading  does  not  appear  to  be  a  highly  sensitive 
function  of  the  manner  in  which  the  wake  spacing  is 
assumed  to  vary  during  the  transient  interval.  Since 
loadings  produced  in  forward  flight  do  not  appear  to  be 
very  sensitive  to  small  variations  in  the  assumed  value 

of  the  vertical  transport  velocity,  it  appears  that  the 
question  of  how  and  when  the  wake  spacing  changes  during 
the  transient  interval  is  not  crucial  to  the  task  of  making 
reasonable  predictions  of  transient  air  loadings.  This 
conclusion  holds  only  for  the  range  of  conditions  studied, 
namely,  maximum  advance  ratio  and  maximum  change  in 
collective  pitch  equal  to  0.  287  and  2.  7  degrees  respectively, 
and  may  not  be  valid  for  more  severe  maneuvers  or  higher 
advance  ratios. 

(7)  Two  deficiencies  of  the  developed  TBL  model  discovered 
in  the  process  of  conducting  this  study  are  the  assumptions 
that  (1)  the  trailing  vortices  trail  free  from  the  one-quarter 
chord  and  (2)  maximum  circulation  on  a  rotor  blade  occurs 
at  a  specific  spanwise  position  on  the  blade.  Calculations 
have  shown  that  the  former  assumption  leads  to  computational 
difficulties  whenever  the  control  point  in  a  given  spanwise 
segment  falls  outside  the  horseshoe  vortex  associated  with 
that  segement.  With  respect  to  the  second  assumption, 
results  computed  for  certain  flight  conditions  demonstrate 
that  the  location  of  maximum  circulation  shifts  markedly 
along  the  span  with  changing  azimuth  and,  accordingly,  the 
program,  as  presently  constituted,  assigns  incorrect 
circulation  strengths  to  the  rolled  up  tip  and  root  vortices. 
These  problems  can  be  eliminated,  respectively,  by  trailing 
vorticity  from  the  trailing  edge  of  the  blade  and  using 
program  logic  that  seeks  out  the  maximum  circulation  for 
assignment  to  the  rolled-up  wake. 
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RECOMMENDATIONS 


The  demonstrated  ability  for  computing  both  periodic  and  nonperiodic 
air  load  distributions  on  rotary  wings  holds  forth  the  promise  of, 
ultimately,  being  able  to  compute  (1)  the  performance,  (2)  the  maneuver 
response,  (3)  the  vibratory  loading  and  stresses,  and  (4)  the  response 
to  gust  disturbances  as  modified  and  controlled  by  the  spiral  wake  that  is 
peculiar  to  the  rotary  wing.  It  should  be  noted  that  this  "promise" 
appears  to  be  more  constrained  by  practical  limitations  on  the  develop¬ 
ment  of  suitable  computing  procedures  rather  than  by  inherent  limitations 
on  the  development  of  a  mathematical  analog  of  the  phenomenon. 
Accordingly,  it  is  recommended  that  continuing  investigations  of  the 
dynamics  of  rotary  wings  pay  particular  attention  to  eliminating 
unnecessary  details  in  the  simulation  while  seeking  to  incorporate  those 
refinements  in  the  theory  essential  to  obtaining  improved  agreement 
between  prediction  and  measurement. 

Although  many  different  approximations  have  been  made  in  the  develop¬ 
ment  of  the  TBL  computational  model,  they  may  be  grouped  in  three 
general  categories:  (1)  aerodynamic  simplifications,  (2)  wake-modeling 
simplifications,  and  (3)  blade-structure  and  rotor -kinematics  simplifica¬ 
tions.  Questions  still  remain  as  to  which  of  these  three  categories 
of  simplifications  constitute  the  primary  source  disagreement  currently 
being  obtained  between  prediction  and  measurement.  Accordingly,  it 
is  recommended,  that  wherever  possible,  improvements  in  the 
computational  model  be  made  sequentially  in  order  to  assess  the  degree 
to  which  a  given  change  in  the  model  influences  its  level  of  accuracy. 

It  is  reemphasized  that  this  recommendation  is  motivated  by  a  strong 
appreciation  for  the  programming  difficulties  and,  more  importantly, 
the  computing  time  and  costs  that  are  associated  with  predictions  of 
SSBL  and  TBL  phenomena. 

The  logical  extension  to  the  analysis  reported  herein  would  be  the 
introduction  of  additional  degrees  of  freedom  such  that  blade  loadings  and 
motions  in  the  flapwise,  pitchwise,  and  edgewise  directions,  together 
with  the  angular  velocity  response  of  the  rotor,  can  be  predicted  as  a 
function  of  transient  inputs  of  blade  pitch  and/or  shaft  torque.  The 
research  objective  recommended  for  this  study  would  be  the  determination 
of  transient  response  characteristics  of  a  rotary  wing  as  influenced  by 
(1)  the  peculiar  blade-wake  relationships  that  exist  or  (2)  the  inertial 
and/or  aerodynamic  coupling  forces  and  moments  that  prevail  during 
the  transient  interval  associated  with  changes  in  blade-pitch  control  or 
changes  in  mechanical  torque  applied  to  the  rotor  shaft.  In  the  latter 
instance,  the  generated  mathematical  model  would  serve  as  a  means 
for  gaining  additional  insight  into  the  dynamics  of  rotors  following  a  power 
failure,  for  example.  Ultimately,  consideration  could  and  should  be 
given  to  including  shaft  and  fuselage  degrees  of  freedom  in  a  TBL 
model  for  purposes  of  investigating  the  extent  to  which  improved 
performance  and  stability  and  control  analyses  can  be  made  as  a  result 
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of  eliminating  the  traditional  assumption  of  uniform  inflow.  Mention 
should  also  be  made  of  the  urgent  requirements  for  performing  research 
that  will  define  the  aerodynamic  behavior  of  rotors  at  very  high  or 
infinite  advance  ratios.  Specifically,  the  high  advance  ratios  to  be 
attained  by  compound  and  stoppable. rotor  helicopter  designs  will 
cause  rotor  blades  to  encounter  flow  conditions  in  which  existing  methods 
of  aerodynamic  analysis  are  no  longer  valid. 
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PRIOR  STUDIES  OF  TRANSIENT  AIR  LOADS  ON  ROTARY  WINGS 


Of  the  total  literature  dealing  with  the  transient  loadings  experienced 
by  rotary  wings,  only  that  portion  in  which  efforts  are  made  to  derive 
predictive  analytical  models  are  of  concern  here.  A  review  has  shown 
that  in  every  case  in  which  analytical  predictions  of  transient  loadings 
have  been  made,  these  predictions  have  dealt  with  the  integrated 
behavior  of  the  rotor  rather  than  with  "microscopic"  detail,  namely, 
the  time  and  3pace  distribution  of  induced  velocity  and  the  resulting 
air  loadings. 

The  efforts  of  Carpenter  and  Fridovich  [13]  and  Rebont,  et  al,  [14] 

(as  mentioned  in  the  Introduction)  are  particularly  worthy  of  mention. 

In  the  first  work,  equations  were  developed  to  predict  the  thrust, 
induced  velocity,  and  flapping  response  of  a  rotor  to  a  rapid  increase 
of  blade  pitch.  A  combination  of  blade-element  theory  and  the  momentum 
theory  was  used,  with  the  momentum  equations  modified  to  apply  to  a 
transient  condition  where  the  inflow  is  being  accelerated.  This 
modification  was  implemented  by  considering  the  mass  of  air  being 
accelerated  to  be  equal  to  the  "apparent  additional  mass"  of  fluid 
associated  with  an  accelerating  impervious  disc.  The  set  of  equations 
yielding  the  desired  representation  of  the  selected  physical  phenomenon 
were  as  follows: 

(1)  Instantaneous  thrust  yielded  by  blade-element  theory 
equated  to  the  instantaneous  thrust  yielded  by  modified 
momentum  theory. 

(2)  The  flapping -motion  equation,  with  blades  assumed 
rigid. 

(3)  The  helicopter  motion  or  lift  equation. 

Experiments  performed  on  a  test  tower  yielded  results  which  were 
compared  with  computed  solutions  of  the  first  two  equations,  the  third 
equation  not  being  needed  since  the  vertical  velocity  of  the  hub  was 
zero  due  to  its  fixed  position.  The  concept  of  apparent  additional 
mass  appeared  to  account  satisfactorily  for  the  measured  average 
inflow  and  thrust  response  to  change  in  collective  pitch. 

The  second  study  [14]  was  similar  to  the  investigation  described 
above.  This  latter  investigation  was  concerned  with  the  response  to  a 
change  in  collective  pitch  in  descending  flight  (near  the  autorotation 
regime)  whereas  the  former  study  was  concerned  with  vertical  takeoff, 
i.  e. ,  with  zero  vertical  velocity  and  zero  thrust  for  the  initial  conditions. 
The  concept  of  apparent  mass  was  similarly  employed  to  modify  the 
momentum  equation  to  yield  instantaneous  values  of  thrust.  Flapping 
was  neglected  and  the  blade-element  thrust  was  combined  with 
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momentum  thrust  to  yield  a  differential  equation  having  the  form  of 
Riccati's  equation.  With  collective  pitch  given  as  a  known  function,  the 
resulting  equation  was  numerically  integrated  for  a  range  of  rotor 
operating  conditions  yielding  results  that  bear  some  resemblance  to 
those  obtained  by  Carpenter  and  Fridovich.  Data  obtained  in  experi¬ 
ments  appeared  to  verify  (qualitatively)  predictions  made  for  various 
nondimensional  rotor  parameters;  however,  these  investigators 
concluded  that  the  apparent  mass  defined  by  Carpenter  and  Fridovich 
should  be  doubled  in  order  to  account  satisfactorily  for  their  experi¬ 
mental  results.  Both  groups  of  investigators  showed  that  the  thrust 
response  is  highly  sensitive  to  the  rate  of  change  of  the  collective 
pitch. 

Michel,  et  al,  [15]  used  the  principle  of  superposition  to  separate  the 
response  of  a  helicopter  to  a  gust  from  the  periodic  response  that 
exists  for  an  articulated  rotor  in  forward  flight.  Equations  and  solution 
techniques  were  developed  to  yield  the  transient  response  of  a  rotor 
upon  entering  a  sharp-edged  gust.  Aerodynamic  loadings  were  computed 
on  the  basis  of  a  time  invariant,  uniform  inflow.  It  should  be  noted  that 
this  study  was  concerned  with  procedures  yielding  overall  gust  load  or 
alleviation  factors  rather  than  with  the  details  of  blade  loading,  blade 
motion,  and  blade  stress.  A  subsequent  investigation  [16]  constituted 
an  experimental  check  on  these  theoretical  results  and  demonstrated 
order  of  magnitude  agreement  between  theory  and  experiment.  The 
experimental  results  indicated,  however,  that  the  blade  flapping 
and  moment  response  is  influenced  by  the  collective  pitch  setting  of 
the  rotor,  in  contrast  to  the  assumption  made  in  the  analysis  [15], 
namely,  that  the  gust  response  is  independent  of  the  steady-state 
periodic  response  that  exists  in  forward  flight. 

Except  for  the  absence  of  means  for  computing  the  detailed  influence 
of  a  vortical  wake,  Gessow  and  Crim  [17]  developed  a  computing 
program  (in  1955)  that  bears  a  close  resemblance  to  the  program  devel¬ 
oped  in  the  investigation  reported  herein.  This  reference  is  cited  last 
since  it  represents  the  state  of  the  art  prior  to  the  present  study  for 
computing  the  transient  flapping  of  helicopter  rotor  blades.  In  fact, 
the  procedure  as  developed  by  Gessow  and  Crim  (or  modified  form) 
is  currently  employed  by  the  helicopter  industry  to  compute  both  the 
steady -state  and  transient  flapping  behavior  of  rigid -blade  rotors.  The 
computational  procedure  is  employed  mainly  for  (1)  investigating 
questions  of  rotor  stability  at  high  advance  ratios,  (2)  determining 
steady-state,  tip -path -plane  orientations,  and  (3)  predicting  perform¬ 
ance  characteristics.  Provision  is  generally  made  for  determining 
blade  thrust  as  a  function  of  lift  and  drag  data  obtained  from  two- 
dimensional  airfoil  tests  over  a  zero-to-360-degree  range  of  angle  of 
attack.  However,  no  provision  is  made  for  introducing  the  influence 
of  temporal  and  spatial  variations  in  induced  velocity,  since  such 
variations  require  that  a  suitable  mathematical  model  of  a  rotor  wake 
be  available. 
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It  should  be  again  pointed  out  that  the  present  study  is  primarily 
concerned  with  the  elimination  of  the  traditional  assumption  of  uniform 
inflow  such  that  the  impulsive  and/or  high  frequency  components  of 
blade  loading  due  to  the  presence  of  the  spiral,  vortical  wake  can  be 
determined  with  reasonable  accuracies.  Previous  research  efforts 
[12]  have  demonstrated  that  this  objective  can  be  achieved  for  SSBL 
calculations,  albeit  at  the  cost  of  considerable  program  complexity 
and  computing  time. 
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REPRESENTATION  OF  THE  ROTOR  WAKE 


GENERAL  REMARKS 


For  purposes  of  computing  the  transient  air  loadings  caused  by  changes 
in  collective  pitch,  the  rotor  and  its  wake  are  modeled  (mathematically) 
in  essentially  the  same  format  that  was  adopted  for  the  SSBL  calcu¬ 
lations  described  in  Reference  12.  Namely,  each  blade  of  the  rotor  is 
represented  by  a  segmented  lifting  line  (bound  vortex)  located  along  the 
quarter  chord  of  the  blade,  and  the  wake  is  represented  by  a  mesh 
(or  lattice  work)  of  segmented,  straight  vortices,  with  the  strength  of 
the  vortex  being  constant  over  the  length  of  the  segment.  Whereas, 
previously,  the  wake  was  considered  to  consist  of  a  large  number  of 
closed  "boxes"  of  vorticity,  constituting  a  convenient  means  for 
"lumping"  the  distributed  vorticity  in  the  wake  into  discrete  elements 
of  shed  and  trailing  vorticity,  the  wake  adopted  for  the  TBL  model  is 
assumed  iO  consist  only  of  trailing  elements.  Omission  of  the  shed 
elements  means  that  we  have  elected  to  ignore  the  influence  of  unsteady 
aerodynamics  (associated  with  the  changing  circulation  on  a  two- 
dimensional  segment  of  the  rotor  blade)  in  computing  air  loadings  in  the 
TBL  calculation.  This  simplifying  assumption  was  adopted  as  a  result 
of  (1)  having  performed  a  study  of  means  for  simplifying  the  wake 
model  (which  study  will  be  described  below),  (2)  recognizing  that 
theoretical  means  for  including  unsteady  aerodynamic  effects  in  the 
TBL  model  are  not  clearly  established,  and  (3)  the  large  motivation 
that  exists  for  reducing  the  complexity  of  the  wake  calculations  when¬ 
ever  it  appears  to  be  reasonably  justifiable.  A  most  important  feature 
of  the  wake  model  that  has  been  employed  heretofore  for  making 
calculations  of  air  loadings  on  rotary  wings  is  the  requirement  to 
specify  beforehand  the  gross  geometry  of  the  vortex  elements  in  the 
wake.  In  practice,  this  has  meant  that  the  steady-state  wake  is 
treated  as  a  "rigid"  wake  in  which  the  pitch  of  the  spiral  is  determined 
by  a  "transport"  velocity  equal  to  the.  mean  velocity  induced  by  the 
rotor,  as  yielded  by  momentum  theory.  A  transient  wake,  equivalent 
to  the  rigid  wake  produced  in  steady,  forward  flight,  is  one  in  which  the 
vorticity  is  trailed  from  a  blade  executing  nonperiodic  rather  than 
periodic  motions  in  space.  Note  that  a  tip-path-plane  is  not  defined. 
However,  the  vorticity,  subsequent  to  being  trailed  by  the  blade, 
acquires  a  geometry  governed  by  a  priori  statements  about  the  steady 
and  time-varying  components  of  the  velocity  transporting  the  vorticity 
in  a  direction  perpendicular  to  the  horizontal  line  of  flight. 

This  a  priori  specification  of  the  time -varying  transport  velocity 
greatly  simplifies  the  model  of  a  rotor -wake  system.  In  so  doing,  we 
eliminate  the  necessity  to  compute  the  self-induced  effects  of  the  wake, 
either  in  terms  of  how  these  self-induced  effects  control  (1)  the  "steady" 
distortion  that  exists  in  a  steady  (periodic)  wake  or  (2)  the  time -varying 
distortion  or  geometry  that  exists  for  a  nonsteady  (nonperiodic)  wake. 
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Thus  we  need  compute  only  the  velocities  induced  by  the  wake  at  the 
rotor  blades,  having  eliminated  the  requirement  for  computing  velocities 
elsewhere  in  the  wake  space. 

Piziali  and  DuWaldt  [10,  11  ]  have  shown  that  a  wake  spacing  (i.  e. ,  of 
the  skewed  spiral  assumed  to  correspond  to  the  wake  of  a  rotary  wing 
in  forward  flight)  based  on  a  downward-transport  velocity  equal  to  the 
momentum -theory  value  of  induced  inflow  yields  an  azimuthal  and  span- 
wise  variation  of  loading  in  reasonable  agreement  with  experiment.  It 
should  be  recognized  that  this  particular  simplifying  assumption  (adopted 
for  the  SSBL  study)  is  the  major  feature  of  the  wake  model  that  results 
in  the  approach  being  a  practicable  one  from  the  standpoint  of  tolerable 
computing  complexity  and  cost.  Accordingly,  no  consideration  was  given 
towards  eliminating  this  assumption  in  this  study  of  transient  *  ladings. 
Note,  however,  that  the  existence  of  a  transient  wake  presents  added 
reasons  for  the  removal  of  arbitrariness  with  respect  to  wake  geometry, 
and  accordingly,  for  the  introduction  of  wake  geometry  as  a  dependent 
variable  in  the  problem.  Approximate  means  for  representing  the 
variable  geometry  in  the  wake,  as  exists  during  the  transient  interval, 
are  discussed  later  in  this  report. 

Although  the  representations  of  wake  vorticity  and  geometry,  as  used  in 
the  SSBL  and  TBL  models,  are  mathematically  equivalent,  there  are 
basic  differences  in  the  solution  format.  In  the  case  of  the  SSBL  calcu¬ 
lation,  an  axis  system  is  placed  in  the  specified  (i.  e. ,  known)  tip-path 
plane.  A  set  of  simultaneous  algebraic  equations  are  derived  in  which 
the  unknown  variables  to  be  determined  are  the  bound  vorticities  that 
exist  on  each  spanwise  Made  segment  at  a  finite  number  of  azimuth 
intervals,  this  number  being  taken  sufficiently  large  to  describe 
adequately  the  azimuthal  variation  of  loading.  (This  procedure  means 
that  the  use  of  10  spanwise  segments  and  24  azimuth  intervals,  for 
example,  requires  the  determination  of  240  unknown  circulation 
strengths.)  in  the  TBL  calculation,  time-varying  loadings  cannot  be 
obtained  from  a  simultaneous  solution  of  equations  that  are  based  on  the 
periodicity  inherent  to  the  steady-state  loading  process.  Rather  the 
bound  vorticity  existing  on  the  blades  at  every  discrete  interval  of  time 
(corresponding  to  an  equivalent  azimuthal  interval)  must  be  determined 
as  a  function  of  the  total  wake  that  exists  at  that  specific  interval  of  time 
Since  at  any  given  time  interval,  the  number  of  unknown  circulation 
strengths  that  must  be  determined  is  the  product  of  the  number  of  span- 
wise  segments  and  the  number  of  blades,  it  is  seen  that,  with  respect 
to  the  requirement  for  solving  a  set  of  simultaneous  equations,  the 
TBL  calculation  is  considerably  less  involved  in  this  regard.  On  the 
other  hand,  whereas  the  SSBL  calculation  requires  that  a  solution  of 
a  very  large  set  of  equations  be  obtained  only  once,  the  mathematical 
analog  of  the  TBL  phenomenon  requires  that  the  equation  set  be  solved 
for  each  time  interval  for  as  many  time  intervals  as  are  required  to 
define  the  transient  response  to  a  given  input  of  blade-pitch  control. 

This  particular  feature  of  the  TBL  calculation  requires  that  the  induced- 
velocity  coefficients  of  the  wake  (as  defined  in  Reference  12)  be 
recomputed  at  each  time  interval.  Note  that  a  SSBL  calculation  requires 
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only  a  single  determination  of  these  coefficients.  Accordingly,  a  TBL 
calculation  procedure  must  pay  a  significant  penalty  in  computing  time 
as  a  consequence  of  the  nonexistence  of  periodicity  in  the  phenomenon. 

The  development  of  the  TBL  calculation  procedure  to  its  present  state 
has  required  the  adoption  of  the  following  simplifying  assumptions  in 
addition  to  the  above-mentioned  (1)  neglect  of  shed  vorticity  and  (2) 
requirement  for  specifying  wake  geometry: 

(3)  The  continuous,  spanwise  distribution  of  bound  vorticity 
existing  on  a  real  blade,  as  well  as  its  continuous 
variation  with  time,  can  be  satisfactorily  defined  by 
using  a  sufficient  number  of  blade  segments  and  azimuth 
intervals,  over  which  segment  and  interval  the  vorticity 
is  assumed  to  be  constant. 

(4)  The  vortices  trailed  into  the  wake  are  equal  to  the 
differences  in  the  strengths  of  the  bound  vorticity  on 
adjacent  segments  of  the  blade. 

(5)  The  effects  of  viscous  dissipation  on  the  wake-vortex 
strengths  can  be  neglected. 

(6)  The  angle  of  attack  of  a  specific  blade  segment  is  the 
"effective"  angle  of  attack  computed  at  the  three-quarter 
chord  point. 

(7)  The  circulation  trailed  into  the  wake  behind  a  stalled 
section  is  the  circulation  corresponding  to  the  lift 
produced  at  the  stall  angle. 

(8)  The  velocities  induced  by  a  rotor  wake  can  be  adequately 
determined  by  application  of  the  Biot-Savart  law. 

(9)  It  is  necessary  to  compute  only  that  component  of 
velocity  induced  perpendicular  to  the  horizontal  line  of 
fligh^the  lateral  and  longitudinal  components  being 
deemed  to  be  relatively  unimportant. 

(10)  It  is  permissible  to  neglect  the  bound  vorticity  on  the 
blades  in  computing  the  velocities  that  are  induced  in 
the  rotating  coordinate  system  since  the  influence  of 
bound  vorticity  is  negligible  at  zero  forward  velocity 
and  remains  negligible  for  low  solidities  and  low  advance 
ratios. 


It  should  be  recalled  that  these  coefficients  are  functions  only  of 
wake  geometry  and  rotor  position. 
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(11)  The  rotor  blades  are  structurally  rigid  and  have  only 
a  flapping  degree  of  freedom. 

(12)  The  rotor  revolves  at  a  constant  angular  velocity 
throughout  the  transient  interval. 

(13)  The  hub  of  the  rotor  continues  to  translate  in  level, 
forward  flight  at  constant  speed  during  the  short 
time  interval  of  interest. 


Further  simplifications  that  have  been  introduced  into  the  analysis 
pertain  to  the  computation  of  the  effective  aerodynamic  angle  of  attack 
and  the  total  aerodynamic  force  created  at  a  given  blade  segment.  For 
example,  the  TBL  calculation  (in  its  present  form)  linearizes  the 
relationships  that  exist  between  circulation  and  induced  velocity  by 
assuming  that  (1)  the  airflow  angle,  <p  ,  is  small,  namely,  that 


(f  *  ta.n~' 


Ur 


. —  Up 

“  ~u7 ' 


where  up  and  Ur  are  as  defined  in  the  Liot  of  Symbols,  and  (2)  lift  is 
linearly  related  to  the  angle  of  attack  up  to  an  assumed  angle  of  stall, 
above  which  the  lift  coefficient  and  circulation  remain  constant.  This 
linearization  was  incorporated  into  the  SSBL  model  described  in 
Reference  12  and  has  been  retained  here  for  purposes  of  simplifying 
and  expediting  the  solution  of  the  equation  set  yielding  the  unknown 
bound  vortices. 


DEVELOPMENT  OF  A  SIMPLIFIED  WAKE  MODEL 

As  part  of  the  overall  effort  to  develop  a  method  for  predicting  the 
nonperiodic  air  loadings  on  a  helicopter  rotor,  considerable  attention 
was  directed  first  towards  achieving  means  for  simplifying  the 
previously  developed  wake  representation.  The  motivation  was  two¬ 
fold.  First,  it  was  essential  that  the  computation  time  associated 

with  computing  the  velocity  coefficients  of  the  wake  be  minimized  - 

particularly  in  view  of  the  additional  computing  that  is  introduced  by 
the  step-by-step  procedure  associated  with  the  TBL  process.  Second, 
it  appeared  desirable  to  introduce  modifications  to  the  original  format 
in  order  to  obtain  a  more  realistic  simulation  of  the  rolling-up  process 
that  occurs  in  the  wake  developed  by  a  rotary  wing.  To  this  end,  an 
investigation  was  made  to  determine  the  extent  to  which  the  original 
wake  representation,  diagrammed  schematically  in  Figure  1,  could  be 
simplified  and/or  modified. 

The  following  wake  modifications  were  examined  for  their  influence 
on  the  induced -velocity  and  loading  distributions  occurring  in  forward 
flight: 
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(1)  Elements  of  shed  vorticity  in  the  wake  were  neglected 

at  progressively  closer  points  to  the  blades  while  retaining 
a  complete  mesh  of  trailing-vortex  segments;  see  Figure 

l  (a). 

(2)  Beyond  the  shed-vortex  truncation  point  of  (1)  above,  the 
trailing-vortex  segments  were  concentrated  ("rolled  up") 
into  a  tip  and  root  vortex;  see  Figure  2  (b).  Provision  was 
also  made  for  positioning  the  root  and  tip  vortices  to 
simulate  wake  contraction. 

(3)  Beyond  the  shed-vortex  truncation  point  of  (1)  above,  the 
trailing-vortex  segments  were  concentrated  into  a  single 
tip  vortex;  see  Figure  2  (c). 

Calculations  were  made  for  the  two-bladed  HU-1A  helicopter  flying  at 
110  knots.  Comparisons  were  made  with  respect  to  the  azimuthal 
variations  in  loading  and  induced  velocity  and  the  harmonic  content  of 
these  variables  up  to  the  tenth  harmonic.  Typical  results  (aerodynamic 
loading  plotted  versus  azimuthal  position  for  a  specific  spanwise 
station)  are  shown  in  Figure  3  for  the  following  wake  configurations: 

(1)  A  complete  grid  of  trailing  and  shed  vorticity. 

(2)  Shed  vorticity  truncated  at  one-quarter  revolution 
behind  blade,  grid  of  trailing  vorticity  retained. 

(3)  No  shed  vorticity,  grid  of  trailing  vorticity  retained. 

(4)  Shed  vorticity  truncated  at  one-quarter  revolution 
behind  blade,  trailing  vorticity  rolled  up  into  a  root 
and  tip  vortex  one-quarter  revolution  behind  the 
blade. 

(5)  Shed  vorticity  truncated  fifteen  degrees  behind 

blade;  trailing  vorticity  "rolled  up"  into  tip  vortex  only. 

Shown  also  in  Figure  3  is  representative  of  similar  plots  made  for 
other  spanwise  stations  and  collectively,  these  plots  demonstrate  that 
the  shed  vorticity,  as  represented  by  the  wake  model  used  in  the  study 
reported  in  Reference  12,  does  not  significantly  influence  the  resulting 
loadings  on  the  blades.  The  difficulties  encountered  in  adequately 
representing  the  shed  vortex  sheet  by  discrete  vortices  located  at  equally 
space  azimuth  intervals  behind  the  blade  are  reviewed  and  discussed 
in  Reference  19. 

Shown  also  in  Figure  3  are  experimental  loadings  obtained  with  the 
HU-1  A  helicopter  [18],  It  appears  that  a  wake  model  consisting  only 
of  trailing-vortex  elements,  trailed  as  a  consequence  of  a  spanwise 


14 


variation  in  bound  vorticity,  does  account  for  a  major  portion  of  the 
loading  phenomena  observed  in  forward  flight.  Since  considerable 
questions  still  remained  as  to  the  adequacy  and  necessity  of  means  for 
representing  the  shed-vortex  sheet,  a  decision  was  to  drop  shed 
vorticity  from  the  wake  model  incorporated  into  the  TBL  program. 

Some  reservation  must  be  retained,  however,  as  to  the  adequacy  of 
this  simplification  when  very  rapid  changes  of  bound  vorticity  are 
produced  as  a  result  of  (1)  very  rapid  changes  in  collective  pitch,  and/ 
or  (2)  very  rapid  changes  in  the  induced-velocity  field.  Admittedly, 
more  calculations  are  in  order  to  assess  the  quantitative  role  of  shed 
vorticity  in  those  flight  conditions  in  which  unsteady  aerodynamic 
effects  would  be  accentuated. 

During  the  wake-simplification  study,  a  few  calculations  were  made  to 
investigate  the  influence  of  wake  contraction  on  the  air  loadings  produced 
on  a  HU-1A  rotor  at  an  advance  ratio  of  0.  18.  These  forward-flight 
calculations  were  inconclusive  in  that  no  definitive  indications  could 
be  seen  as  to  whether  the  simulation  of  the  contraction  in  the  rolled-up 
wake  yields  improved  correlation  with  experiment.  Subsequent 
calculations  made  for  hovering  flight,  on  the  other  hand,  produced 
just  the  opposite  result,  namely,  that  simulation  of  wake  contraction 
is  a  feature  of  a  wake  model  that  is  essential  to  the  prediction  of  the 
aerodynamic  loading  in  hover.  A  few  remarks  relative  to  this  point  are 
in  order. 

Since  the  objective  of  the  SSBL  research  program  was  the  development 
of  means  for  predicting  the  highly  variable  ioadings  produced  in  forward 
flight,  no  attention  had  been  directed  to  hovering  flight.  Thus,  no 
knowledge  was  available  as  to  the  ability  of  the  SSBL  wake  model  to 
predict  total  loads  and  load  distributions  of  reasonable  accuracy  for  a 
hovering  rotor.  In  this  investigation  the  response  of  a  hovering  rotor 
to  rapid  changes  of  collective  pitch,  was  however,  a  matter  of  interest. 

It  was,  therefore,  necessary  to  obtain  a  steady-state  loading  in  hover 
for  purposes  of  representing  the  initial  condition  that  prevails  prior  to 
the  initiation  of  the  transient  input.  Calculations  showed  that  a  wake 
model  without  any  contraction  assumed  for  the  rolled-up  portion  yielded 
both  a  total  load  and  spanwise  distribution  of  loading  that  were  in  poor 
agreement  with  the  results  measured  on  an  H-34  helicopter  [20]. 

Figure  4  shows  the  nature  of  this  discrepancy  and  presents  also  the 
considerable  improvement  in  the  calculated  result  when  the  rolled-up 
tip  vortex  is  assumed  to  be  located  at  r/R  =  0.  825  and  the  rolled-up 
root  vortex  is  assumed  to  be  located  at  r/R  -  0.425. 

It  is  concluded  that  realistic  predictions  of  loading  during  hover  can  be 
obtained  only  if  the  wake  geometry  is  suitably  specified  to  approximate 
the  actual  contraction  of  the  wake.  Further,  the  limited  wake  studies 
conducted  in  this  investigation  indicate  that  the  variable  loadings 
produced  for  forward  flight  in  the  vicinity  of  fji-  0.  18  are  caused  by 
geometrical  features  of  the  wake  that  can  be  described  satisfactorily 
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irrespective  of  small  errors  in  the  estimates  of  (1)  the  degree  of  wake 
contraction  and  (2)  the  value  of  the  vertical-transport  velocity.  The 
spanwise  load  distributions  produced  during  hovering,  on  the  other 
hand,  are  sensitive  functions  of  geometry  of  the  wake  —  in  particular, 
the  pitch  of  the  vortex  spiral  and  the  diameter  of  the  spiral  formed  by 
the  rolled-up  tip  vortex. 

THE  COMPUTATIONAL  MODEL 

A  mathematical  analog  of  an  articulated  rotor,  with  rigid  blades  flapping 
in  the  induced-velocity  field  created  by  the  vortical  wake  of  the  rotor, 
consists  of  several  major  parts.  These  are  as  follows: 

(1)  The  equations  that  yield  the  geometry  of  the  wake,  i.  e. , 
the  spatial  location  and  orientation  of  elements  of  vorticity 
in  the  wake  relative  to  the  blades  of  the  rotor  at  each  and 
every  instant  of  time.  (Wake  geometry  is  determined  by 
the  specified  flight  condition  and  rotor  configuration,  the 
motions  of  the  rotor  blades,  and  the  vertical. transport 
velocity.  This  latter  variable  is  specified  as  a  constant 
for  calculating  steady-state  loadings  am  as  a  function  of 
time  when  computing  transient  loadings. ) 

(2)  The  application  of  the  Biot-Savart  law  to  compute  the 
induced-velocity  coefficients  corresponding  to  the  wake 
geometry  existing  at  a  given  instant  of  time.  (These 
induced-velocity  coefficients  represent  the  velocity 
induced  by  a  given  segment  of  wake  at  a  specific  location 
in  the  rotor  disc  per  unit  strength  of  vorticity  in  the  wake 
element.  ) 

(3)  The  assignment  of  circulation  strengths  to  each  element 
of  wake  at  each  instant  of  time  as  a  function  of  all  the 
vorticity  that  has  been  trailed  at  earlier  instants.  These 
vorticities  are  functions  of  the  bound  vorticities  that 
existed  at  earlier  times. 

(4)  The  summing  of  induced-velocity  contributions  from  all 
elements  of  the  wake  except  for  the  elements  of  unknown 
strength  immediately  adjacent  to  the  blades. 

(5)  The  equations  defining  the  angle  of  attack  of  the  blade 
segments  as  a  function  of  the  kinematics  of  the  rotor  and 
the  velocities  induced  by  all  of  the  vorticity  in  the  wake, 
except  the  vorticity  being  trailed  at  the  given  instant  of 
time. 

(6)  The  set  of  simultaneous  equations  that  relate  at  the  given 
instant  of  time  the  unknown  bound  circulations  (and 
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likewise  the  trailing  vorticities  adjacent  to  the  blades) 
to  those  values  of  circulation  corresponding  to  the 
angles  of  attack  determined  in  (5)  above. 

(7)  The  equations  defining  aerodynamic  loading,  moment  about 
the  flapping  hinge,  and  angular  acceleration  of  the  blade 
about  the  flapping  hinge. 

(8)  A  numerical  integration  procedure  to  obtain  the  angular 
velocity  and  displacement  of  the  blade  about  the  flapping 
hinge. 

Figure  5  is  a  schematic  of  the  above -outlined  mathematical  analog. 

(The  details  of  the  complete  mathematical  model  are  presented  in 
Appendix  1.)  It  is  sufficient  to  note  here  that  the  computations 
indicated  in  Figure  5  are  carried  out  for  as  many  time  intervals  as 
may  be  necessary  to  (1)  establish  a  periodic  loading  and  wake  and  (Z) 
determine  a  complete  time  history  of  the  transient  response  to  a 
specified  change  in  collective  pitch,  that  is,  until  the  loading  and  wake 
again  become  periodic  with  time. 

The  calculations  diagrammed  in  Figure  5  require  that  an  axis  system 
be  defined  fox*  describing  both  the  flapping  response  of  the  rotor  and  the 
geometry  of  the  wake.  Figure  6  shows  the  axis  system  adopted 

for  this  purpose.  Note  that  the  origin  is  placed  at  the  hub  of  an 
articulated  (flapping)  rotor,  with  the  jt  axis  directed  rearward,  parallel 
to  the  horizontal  direction  of  flight.  It  should  be  recalled  that  the 
analysis  is  restricted  to  level-flight  operations  or  to  the  situation  that 
prevails  when  the  rotor  hub  is  held  fixed  in  a  wind  tunnel.  Figure  6 
shows  that  the  rotor  shaft  has  been  assumed  to  be  tilted  only  in  the  x>-y 
plane.  Thus  we  ignore  any  lateral  tilt  of  the  rotor  shaft,  since  this 
tilt  angle  is  usually  small  during  trimmed,  level  flight.  A  conventional 
shaft  axis  system  is  used  to  describe  the  flapping  and  feathering  of  the 
blades  and  the  motions  and  displacements  of  the  blades  (as  seen  in  the 
shaft  axis  system)  are  transformed  to  equivalent  quantities  in  the 
system.  Only  the  j-axis  component  of  the  total  induced-velocity  vector 
is  included  in  the  analytical  model,  as  was  indicated  earlier.  In 
accordance  with  the  adopted  axis  system,  velocities  induced  in  the 
negative  ^-direction  are  designated  as  negative  quantities. 


17 


COMPARISON  OF  THEORY  AND  EXPERIMENT 


STEADY -STATE  LOADINGS 

The  measured  air  loads  used  to  test  the  adequacy  of  the  TBL 
computational  model  for  predicting  periodic  loadings  in  steady  flight 
were  obtained  from  two  sources:  full-scale  flight  tests  of  the  H-34 
helicopter  [21  ]  and  wind-tunnel  tests  of  a  full-scale  H-34  rotor  [22]. 

The  first  set  of  data  is  a  portion  of  the  data  originally  used  to  evaluate 
the  adequacy  cf  the  SSBL  calculation  procedure  [12].  Data  obtained  at 
one  specific  flight  condition  (  /*-  =  0.  18;  V  -  118  fps)  are  employed  herein 
to  examine  the  extent  to  which  the  current  SSBL  and  the  TBL  computa¬ 
tional  procedures  yield  results  in  agreement  with  each  other  and  with 
experiment.  Although  it  is  patently  true  that  the  successful  prediction 
of  transient  loadings  requires,  a  priori,  a  model  that  correctly  predicts 
periodic  loadings  for  steady  flight,  it  was  felt  that  steady -state 
comparisons  are  also  mandatory  to  demonstrate  that  the  TBL  procedure 
yields  periodic  loadings  for  steady  forward  flight  equivalent  to  that 
calculated  by  another  computational  method. 

Figure  7  presents  the  aerodynamic  loadings  at  seven  spanwise  stations 
as  a  function  of  azimuth,  as  obtained  by  measurement  and  by  computa¬ 
tion  with  the  SSBL  and  TBL  programs.  The  differences  between  the 
results  produced  by  the  two  computational  procedures  can  be  attributed 
to  the  differences  that  exist  between  first-harmonic  flapping,  that  is, 
the  tip-path-plane  orientation.  It  should  be  noted  that  the  SSBL  calcu¬ 
lation  procedure  described  in  Reference  12  requires  that  the  tip-path- 
plane  orientation  be  known  and  specified  as  an  input  to  the  program.  In 
the  TBL  calculation,  however,  the  blades  are  free  to  flap  and  the  steady- 
state  orientation  of  the  tip-path  plane  is  a  computed  quantity.  Small 
differences  between  specified  and  computed  first-harmonic  flapping  of 
the  blade  result  in  significant  differences  in  the  first-harmonic 
variation  of  angle  of  attack.  These  differences  in  computed  angle  of 
attack  caused  by  small  discrepancies  in  first-harmonic  flapping  are 
responsible  for  the  discrepancies  in  computed  loadings  seen  in  Figure  7. 
Notwithstanding  this  failure  of  the  SSBL  and  TBL  calculation  procedures 
to  agree  on  first-harmonic  components  of  the  air  loading,  a  comparison 
of  the  measured  and  computed  loadings  in  Figure  7  shows  that  theory 
predicts  the  higher  frequency  variations  in  measured  load  reasonably 
well. 

Figure  8,  a  plot  of  induced  velocities  versus  azimuth  for  seven  different 
spanwise  positions,  demonstrates  that  the  SSBL  and  TBL  calculation 
procedures  yield,  for  all  practical  purposes,  the  same  induced -velocities 
in  the  plane  of  the  rotor  disc.  This  figure  is  included  for  purposes  of 
demonstrati  ig  that  the  differences  between  the  calculated  loadings 
presented  in  Figure  7  cannot  be  attributed  to  a  lack  of  equivalence  in  the 
induced  velocity  field  yielded  by  two  different  calculation  procedures. 
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The  question  as  to  why  the  TBL  computation  yields  a  tip-path-piane 
angle  of  4.2  degrees  in  contrast  to  the  angle  of  3.0  degrees  specified 
as  the  input  to  the  SSBL  computation  cannot  be  answered  in  a  precise 
manner.  There  are  many  approximations  and  simplifications  (as  well 
as  inaccuracies  of  measurement)  that  lead  to  small  differences  in  load 
and  load  distribution  which,  in  the  aggregate,  can  cause  differences  in 
blade  motion  of  the  magnitude  obtained.  It  should  be  emphasized  that 
while  accurate  prediction  of  blade  motions  is  necessary  to  obtain 
accurate  predictions  of  blade  loadings,  it  is  by  no  means  sufficient. 

Figure  9  is  a  comparison  of  calculated  periodic  air  loadings  (i.  e.  ,  in 
steady  forward  flight  with  data  obtained  in  full-ucale  wind-tunnel  tests 
[22],  The  loadings  plotted  in  this  figure  correspond  to  the  initial  steady 
state  achieved  in  Run  1  (see  table  on  the  following  page  for  listing  of 
test  conditions  -  Runs  1,  2,  and  3)  prior  to  the  initiation  of  a  collective- 
pitch  change.  In  general,  the  agreement  between  theory  and  wind-tunnel 
test  appears  to  be  as  good  or  better  than  that  achieved  with  full-scale 
flight  data  (see  Figure  7).  Since  the  agreement  that  is  obtained  is,  in 
part,  a  function  of  the  advance  ratio  and  total  force,  it  is  not  possible 
to  draw  any  firm  conclusions  as  to  whether  the  improved  agreement 
reflects  more  accurate  input  data  or  is  the  result  of  the  air  loadings 
being  less  sensitive  to  the  wake  fo**  this  particular  flight  condition.  It 
should  be  noted,  however,  that  the  comparisons  made  in  Figure  7  involve 
calculations  in  which  no  wake  contraction  was  introduced  whereas  the 
tip  and  root  vortices  were  assumed  to  roll  up  45  degrees  of  azimuth 
behind  the  blade  at  the  0.  90  and  0.  375  radial  stations,  respectively, 
in  the  calculations  plotted  in  Figure  9.  No  calculations  have  been  made 
without  wake  contraction  being  assumed  for  the  test  conditions 
corresponding  to  the  results  presented  in  Figure  9. 


RESPONSE  TO  CHANGE  IN  COLLECTIVE  PITCH 

The  measured  air  loads  used  to  test  the  adequacy  of  the  TBL  computa¬ 
tional  model  for  predicting  the  transient  response  to  a  change  in 
collective  pitch  were  obtained  in  wind-tunnel  tests,  as  indicated  above. 
Data  were  obtained  at  two  "flight"  conditions.  For  the  first  condition, 
one  transient  response  was  produced  by  a  rather  slow  increase  in 
collective  pitch.  In  the  second  condition,  two  transients  were  produced, 
one  in  response  to  a  more  rapid  increase  in  collective  pitch,  a 
second  in  response  to  a  rapid  decrease.  Data  pertaining  to  the  wind- 
tunnel  test  conditions  are  tabulated  in  the  table. * 


Originally  a  more  exhaustive  test  program  was  planned,  but 
instrumentation  difficulties  forced  the  curtailment  of  the 
transient-response  tests. 
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In  order  to  facilitate  a  comparison  between  theory  and  experiment,  it 
was  necessary  first  to  reduce  the  test  data  for  purposes  of  determining 
the  cyclic -pitch  settings  and  time  history  of  collective  pitch.  Since 
blade  pitch  was  measured  on  blades  separated  by  180  degrees  of 
azimuth,  collective  pitch  could  be  resolved  by  summing  the  output  of 
these  pitch  transducers.  Cyclic  pitch  was  then  determined  by 
subtracting  collective  from  the  pitch  of  the  instrumented  blade. 


TABLE  I 

MIND  TUNNEL  TEST  CONDITIONS 


RUN  1 

RUN  2 

RUN  3 

AVERAGE  TEST  VELOCITY,  FT. /SEC. 

131.0 

IB6.0 

IB6.5 

ROTOR  ANGULAR  VELOCITY,  RAD./SEC. 

23.2 

23.2 

23.2 

ADVANCE  RATIO 

0.202 

0.286 

0.287 

SHAFT  ANGLE  OF  ATTACK,  DEG. 

10 

5 

5 

AVERAGE  AIR  DENSITY,  SLUQS/CU.FT. 

0. 002?r  ' 

0.002276 

0.002272 

INITIAL  TOTAL  LIFT,  LBS. 

6710 

8945 

8715 

FINAL  TOTAL  LIFT,  LBS. 

12162 

• 

2037 

TOTAL  COLLECTIVE-FITCH  INCREMENT,  OEQ. 

1.775 

1.607 

-2.663 

RATE  OF  CHANGE  OF  COLLECTIVE  FITCH 
(NOMINAL  VALUE)  DEG./SEC 

6.55 

12.95 

-47.2 

*N0  WIND-TUNNEL  BALANCE  READING  BECAUSE  OF  EXCESSIVE  V I  BRAT  I  ON . 


The  above  procedure  worked  very  well  with  respect  to  determining  the 
relative  variation  in  collective-pitch  settings,  but  the  absolute  values 
of  collective  pitch,  as  indicated  by  the  measured  traces  and  the  supplied 
calibration  data,  were  so  small  that  negative  thrust  would  result  in  a 
subsequent  computation.  It  proved  expedient  to  resort  to  the  charts 
presented  in  Reference  23  to  determine  appropriate  settings  of 
collective  pitch  and  to  utilize  further  that  portion  of  the  TBL  program 
which  operates  with  a  specified  inflow  to  check  these  settings  for 
suitable  accuracy.  Subsequent  to  this  step,  as  made  necessary  by  the 
invalid  experimental  data,  three  complete  transient  responses  were 
computed  for  comparison  with  actual  measurement. 

Figures  10,  11,  and  12  present  comparisons  of  measure  and  predicted 
transient  air  loadings  for  seven  span  positions  for  test  runs  1,  2,  and 
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3,  respectively.  '  Shown  also  are  predicted  and  measured  flapping,  the 
time  history  of  collective  pitch  (as  measured  and  as  used  for  the  input 
to  the  calculation),  and  the  assumed  variation  in  the  mean  induced 
velocity,  W,  governing  the  spacing  of  the  wake  during  the  steady  and 
transient  interval.  The  assumed  timewise  variation  in  w  was  based  on 
the  results  of  hovering  calculations  that  are  discussed  in  the  following 
section  of  the  report.  As  a  result  of  this  hovering  study,  there  is 
reason  to  believe  that  the  computed  air  loadings  would  not  be  modified 
significantly  on  assuming  a  different  time  variation  in  j?.  Instead,  it 
appears  that  wake  distortion  arising  from  self-induced  effects,  in  all 
likelihood,  contributes  errors  in  the  description  of  wake  geometry  that 
are  considerably  more  significant  than  the  errors  that  result  from  the 
very  approximate  means  used  to  simulate  the  time-varying  spacing  of 
the  wake  spiral  during  the  transient  interval. 

Figures  10,  11,  and  12  show  that  the  TBL  model  developed  in  this 
study  simulates  in  remarkably  good  fashion  the  overall  behavior  of  a 
rotary  wing  subjected  to  a  change  in  collective  pitch.  Although  the 
computations  do  not  produce  precise  agreement  with  the  measured 
time  histories  of  air  loadings,  there  are  many  areas  of  qualitative 
and  quantitative  agreement.  It  is  seen  that  certain  portions  of  the 
high-frequency  content  in  the  loadings  are  predicted  by  theory 
remarkably  well.  Note  should  be  taken  of  the  agreement  between  theory 
and  experiment  as  reflected  by  the  different  patterns  of  loading  that 
occur  at  the  various  radial  positions  on  the  blade.  Note  also  the 
agreement  achieved  between  measured  and  predicted  flapping. 

These  comparisons  of  measured  and  calculated  responses  would 
appear  to  indicate  that  continued  refinements  of  the  analysis  with 
respect  to  (1)  the  aerodynamics,  (2)  improvements  in  the  description 
of  wake  geometry,  and  (3)  expanded  modeling  of  the  rotor  structure 
will  lead  to  an  even  better  prediction  capability.  It  is  concluded  that 
the  wake  model,  in  its  present  form,  together  with  a  model  of  the 
kinematics  and  dynamics  of  a  simple  articulated,  rigid-blade  rotor 
simulates  the  essential  components  of  the  TBL  phenomenon  resulting 
from  typical  changes  in  collective  pitch.  Admittedly,  additional  data 
are  needed  at  lower  and  higher  advance  ratios  in  order  to  further 
confirm  or  qualify  this  conclusion. 


It  should  be  noted  that  a  r/R  =  0.  95,  the  calculated  response 
corresponds  to  r//f  =  0.  96  as  a  result  of  a  requirement  (to  be 
discussed  later)  to  compute  the  loadings  at  eight  spanwise 
segments. 
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A  comment  is  in  order  with  respect  to  the  total  thrusts  that  are 
predicted  by  a  nonuniform -inflow  calculation.  As  was  mentioned  above, 
the  collective  pitch  was  adjusted  until  a  uniform-inflow  calculation 
yielded  a  calculated  total  thrust  equal  to  the  total  thrust  measured  in 
the  wind  tunnel.  In  actual  practice,  allowance  was  made  for  the  fact 
that  the  integration  (trapezoidal)  of  load  along  the  span  yields  calculated 
total  blade  loads  that  are  less  than  the  true  load.  For  example,  the  use 
of  eight  spanwise  segments  in  the  calculation  results  in  a  th  ee-percent 
reduction  from  the  load  yielded  by  the  use  of  nine  segments.  Accordingly, 
the  initial  collective  pitch  was  set  in  Run  1  to  produce  an  average  load 
of  8450  pounds,  with  the  uniform-inflow  computation  indicating  a  four- 
per-revolution  oscillatory  component  of  approximately  60  pounds.  It 
was  found  that  the  steady  total  load,  as  yielded  by  a  nonuniform -inflow 
calculation,  is  approximately  8150  pounds  and  has  a  four -per -revolution 
oscillatory  component  of  approximately  755  pounds.  Two  observations 
can  be  made: 

(1)  The  average  total  lift  resulting  from  a  nonuniform -inflow 
calculation  is  less  than  that  yielded  by  a  uniform -inflow 
calculation. 

(2)  A  significant  level  of  oscillatory  total  shear  at  the  rotor 
hub  is  predicted  by  a  nonuniform -inflow  calculation 
whereas  very  little  oscillatory  component  is  predicted 
by  a  uniform -inflow  calculation. 

The  above  comparisons  have  significance  with  respect  to  calculations 
of  rotor  performance  and  transmitted  shears. 

Note  that  Figures  10,  11,  and  12  do  not  readily  indicate  whether 
discrepancies,  if  any,  exist  between  the  calculated  and  measured 
mean  loadings.  The  calculations,  as  performed,  yielded  mean  total 
thrust  approximately  4  percent  less  than  the  initial  and  final  mean 
steady-state  thrusts  measured  for  all  three  runs  with  one  exception  — 
the  calculated  average  thrust  for  the  final  steady  state  achieved  in 
run  3  is  approximately  2600  pounds  (i.  e.  ,  about  30  percent  greater 
than  the  measured  value).  These  discrepancies  in  the  mean  thrust 
levels  obviously  do  not  detract  from  the  basic  validity  of  Figures  10, 

11,  and  12  for  demonstrating  the  ability  of  the  TBL  model  to  predict 
variations  in  air  loading,  either  with  respect  to  spanwise  position  on 
the  blade  or  azimuthal  location  in  the  rotor  disc. 
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ADDITIONAL  NUMERICAL  INVESTIGATIONS 


The  development  of  a  method  for  computing  the  transient  air  loadings 
on  rotary  wings  has  necessarily  resulted  in  several  additional  numerical 
studies.  The  results  obtained  in  these  studies  have  yielded  additional 
insight  on  (1)  the  nature  of  the  blade-loading  phenomenon  and  (2)  the 
capabilities  and  shortcomings  of  the  developed  analytical  procedure. 
Several  of  these  results  will  be  discussed  in  this  section  of  the  report. 

Consideration  of  the  blade -load  prediction  task  shows  that  there  are 
several  sources  of  error.  One  source  is  the  erroneous  contribution  to 
the  angle  of  attack  of  a  blade  segment  that  results  from  geometric 
angles  of  pitch  and  blade  pitching/plunging  velocities  that  are  not  in 
accord  with  actuality  as  a  result  of  assuming  that  the  blades  are  rigid 
or  as  a  result  of  insufficient  rigor  in  modeling  either  the  elastic  structure 
of  the  blade  or  the  kinematics  of  the  hub.  A  second  source  is  the  erron¬ 
eous  induced-velocity  field  that  results  from  our  inability  to  model  the 
wake,  both  its  geometry  and  strength,  with  sufficient  exactitude.  With 
respect  to  the  first  source  of  error,  studies  are  underway  in  several 
organizations  (including  CAL)  to  determine  the  extent  to  which  pre¬ 
dictions  of  blade  loading  can  be  upgraded  by  appropriate  modeling  of  the 
elastic  properties  of  rotor  blades  [19].  With  respect  to  the  second 
source  of  error,  the  extent  to  which  improved  modeling  of  the  geometry 
of  the  wake  will  yield  significant  gains  in  load-prediction  accuracies  is 
a  matter  of  much  concern.  On  the  other  hand,  there  is  considerable 
motivation  (as  was  pointed  out  earlier)  for  simplifying  the  model  of  the 
wake  in  order  to  reduce  the  program  size  and  the  amount  of  computing 
time  that  is  currently  required.  Further,  questions  of  the  adequacy  of 
wake  models  become  considerably  more  involved  when  transient 
(nonperiodic)  loadings  are  of  interest  rather  than  the  periodic  loadings 
produced  in  steady,  forward  flight.  Accordingly,  consideration  is 
given,  in  this  section  of  the  report,  to  certain  aspects  of  wake  modeling, 
both  for  the  steady-state  and  transient  flight  regimes. 

The  contribution  of  the  wake  in  causing  blade-load  variations  of  high 
harmonic  order  is  graphically  demonstrated  in  Figure  13,  Computed 
total  lifts  on  an  H-34  rotor  blade  are  shown  plotted  versus  azimuth.  The 
dashed  curve  is  the  total  load  that  is  computed  when  the  inflow  is  assumed 
to  be  constant.  The  solid  curve  is  the  result  that  is  obtained  when  the 
inflow  is  a  function  of  the  velocities  induced  by  the  vorticity  in  the  wake. 

Although  the  vorticity  distribution  in  the  wake  corresponding  to  the  lift 
on  the  blade  represented  by  the  solid  line  is  significantly  different  from 
that  associated  with  the  lift  on  the  blade  represented  by  the  dashed  line, 
we  find  that  the  induced-velocity  distributions  produced  are  very  similar. 
Figure  14  shows  the  spanwise  variation  of  induced  velocity  on  four  blades 
if  the  computation  is  made  at  an  instant  of  time  at  which  all  the  vortex 
elements  in  the  wake  have  strengths  corresponding  to  the  load  variation 
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represented  by  the  dashed  line  in  Figure  13.  Figure  14  also  shows  the 
induced  velocity  versus  span  that  results  from  a  wake  having  the  vortex 
strengths  associated  with  the  solid  line  in  Figure  13.  It  appears 
reasonable  to  conclude  (at  least  on  a  tentative  basis)  that  the  velocities 
induced  in  the  plane  of  the  rotor  disc  are  more  sensitive  functions  of 
wake  geometry  than  of  the  vortex  strengths.  These  calculations  appear 
to  lend  support  to  the  hypothesis  that  a  wake  having  a  constant  strength 
corresponding  to  the  mean  total  load  on  a  rotor  would  produce  a  good 
approximation  of  the  velocity  field  that  is  induced  in  the  plane  of  the 
rotor  disc.  Additional  numerical  studies  are  required  to  confirm  this 
hypothesis  fully  and  to  establish  the  degree  of  approximation  involved. 

The  problem  of  modeling  rotor -wake  geometry  during  the  transient 
interval  that  is  associated  with  a  change  in  collective  pitch  has  been 
mentioned  previously.  Note  that  it  is  theoretically  possible  to  account 
for  both  the  steady- state  and  transient  distortion  in  a  wake  by  computing 
the  velocities  that  are  induced  everywhere  in  the  wake  space.  However, 
the  realities  of  required  computing  time  and  cost  rule  out  this  approach 
to  the  problem.  Fortunately,  reasonably  good  predictions  can  be  made 
of  the  periodic  loadings  in  steady,  forward  flight  when  the  self-induced 
distortions  in  the  steady  wake  are  ignored  [  1 2  J .  On  the  other  hand,  it 
is  not  possible  to  ignore,  a  priori,  the  time -varying  geometry  of  the 
transient  wake  since  the  assumption  made  with  regard  to  a  uniform  and 
constant  velocity  transporting  wake  elements  in  the  downward  direction 
may  no  longer  be  admissible.  A  short  numerical  investigation  of  this 
point  was  made. 

Since  the  mean  velocity  induced  by  the  wake  is  greatest  in  hovering 
flight,  it  was  hypothesized  that  the  transient  loadings  resulting  from  a 
collective -pitch  change  introduced  during  hover  would  be  most  sensitive 
to  transient  spacing  in  the  wake  spiral.  Accordingly,  a  series  of 
calculations  were  made  for  a  H-34  rotor  at  zero  flight  velocity  in  order 
to  investigate  the  importance  of  transient  wake  effects.  Responses  to 
both  positive  and  negative  changes  in  collective  pitch  were  obtained. 

Figure  15  is  a  time  history  of  loading,  total  thrust  on  rotor,  total 
thrust  on  blade,  flap  angle,  circulation,  and  induced  velocities  as 
caused  by  a  change  in  collective  pitch  of  approximately  four  degrees. 

The  solid  curves  were  obtained  for  a  wake  whose  spacing  remains  fixed 
at  its  initial  steady-state  value.  The  dashed  curves  were  obtained  when 
the  mean  value  of  induced  velocity,  transporting  the  wake,  was  assumed 
to  vary  with  time  as  shown.  *  Note  that  the  collective -pitch  ramp  is 


As  a  result  of  an  error  in  the  computation  of  the  aerodynamic 
moment  about  the  flapping  hinge,  all  hovering  computations  predict 
a  coning  angle  that  is  much  too  large.  Conclusions  drawn  with 
respect  to  loading  and  the  effects  of  transient  wake  geometry  are 
still  valid,  however. 


24 


completed  in  a  time  interval  required  for  the  rotor  to  rotate  through 
60  degrees  of  azimuth.  The  blade-pitch  rates  that  exist  during  the 
input  interval  are  quite  large  and  the  observed  discontinuities  in  the 
loading  time  histories,  etc.  ,  result  from  the  effective  angles  of  attack 
caused  by  the  high  rate  of  blade  pitch  that  exists  over  the  input  interval. 
We  observe  that  the  total  thrust  on  the  rotor  reaches  its  peak  value 
approximately  three-quarters  of  a  revolution  after  initiation  of  the 
collective-pitch  ramp.  It  should  be  noted  that  the  overall  behavior  of 
the  rotor,  i.  e. ,  the  predicted  total  thrust  and  flap  response,  is  similar 
to  the  experimental  results  obtained  previously  by  Carpenter  and 
Fridovich  [13], 

Of  particular  interest  to  this  investigation  are  the  changes  in  loading 
and  load  distribution  that  result  from  a  transient  input  and  subsequent 
variations  in  the  geometry  of  the  wake.  Note  that  the  loadings  do  not 
reach  their  peak  values  simultaneously;  specifically,  the  center-span 
blade  segments  reach  their  peak  loading  prior  to  segments  near  the  tip 
of  the  blade.  * **  This  result  is  caused  by  differences  in  the  timewise 
variation  in  induced  velocity  at  the  various  radial  positions.  The  time 
histories  show  that  the  induced  velocities  change  immediately  following 
the  initiation  of  the  input  at  those  stations  located  at  the  tip  and  root  ends 
of  the  blade.  At  the  median  stations,  there  is  a  time  delay  before  any 
significant  change  is  exhibited  by  the  induced  velocity. 

When  the  mean  induced  velocity,  W ,  varies  as  shown  by  the  dashed  line  , 
we  find,  first,  that  the  final  steady  state  is  attained  at  an  earlier  time 
since  the  increased  wake  spacing  decreases  the  influence  of  the  "older" 
wake  elements.  Second,  we  observe,  as  expected,  that  the  final  total 
thrust  on  the  rotor  is  increased.  Third,  the  time  histories  of  blade 
loading  are  similar  to  those  computed  previously,  with  the  exception 
that  at  r/R  -  0.  55,  0.  65,  and  0.  75  we  observe  (1)  an  "undershoot"  as 
well  as  an  "overshoot"  in  loading,  and  |2)  greater  percentage  differences 
in  the  final  steady- state  loading.  Fourth,  these  observed  changes  in 
load  distribution  along  the  span  are  accounted  for  by  the  different 
induced-velocity  field  that  results  when  the  wake  spacing  is  assumed  to 
increase  in  the  manner  shown.  It  is  seen  that  the  increased  wake 
spacing  causes  a  relatively  small  increase  in  the  velocities  induced 
near  the  tip  of  the  blade,  whereas  at  r/R  -  0.  55,  0.  65,  and  0.  75  a 
significant  reduction  occurs  in  the  velocities  induced  by  the  wake.  In 
general,  it  may  be  stated  that  it  is  primarily  the  mid-span  locations  on 
the  blade  that  exhibit  a  change  in  steady-state  loading  due  to  a  specified 
change  in  wake  geometry  and  these  locations  likewise  exhibit  the  largest 
transient  effects  that  depend  on  the  manner  in  which  the  mean  value  of 
the  induced  velocity  is  assumed  to  change. 


*  The  loading  at  the  25-percent  span  position  remains  constant  since 
this  segment  of  the  blade  is  stalled. 

**  The  final  value  of  7  is  based  on  a  prior  estimate  of  the  final  total 
thrust. 
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The  response  of  a  hovering  rotor  to  a  reduction  in  collective  pitch  is 
presented  in  Figure  16.  The  mean  induced  velocity, py  ,  was  assumed 
to  vary  in  the  manner  plotted.  This  assumption  was  based  on  a  prior 
estimate  of  the  final  value  of  total  thrust  to  be  achieved  with  the  reduced 
collective  pitch,  with  the  timewise  variation  being  based  on  the  induced- 
velocity  responses  computed  for  the  increase  in  collective  pitch. 

Another  aspect  of  the  wake  model,  deserving  of  consideration,  is  the 
assumption  that  the  trailing  vorticity  streams  rearward  (with  the  local 
velocity  vector)  from  the  bound  vortex  located  at  the  one -quarter  chord. 
This  assumption  has  created  a  computational  limitation  in  that  an  upper 
limit  exists  on  the  advance  ratios  that  can  be  treated,  this  limit  being 
a  function  of  the  spanwidths  specified  for  the  blade  segments  used  in 
the  computation.  The  limitation  arises  from  the  "lumping"  process 
used  to  represent  the  trailing  vortex  sheet.  Specifically,  the 
representation  of  the  trailing  wake  by  discrete  vortices  implied 
selection  of  control  points  for  computing  induced  velocities  that 
are  midway  between  trailing  vortex  filaments.  Under  the  present 
assumption  of  vortices  trailing  free  from  the  one-quarter  chord,  it  is 
possible  for  the  three-quarter  chord  point  on  the  blade  to  fall  outside 
the  horse  shoe  vortex  associated  with  that  particular  blade  segment. 

This  resultant  geometry  can  occur  for  a  combination  of  very  narrow 
segments  (usually  taken  at  the  blade  tip)  and  high  advance  ratios.  The 
computation  will  fail  whenever  this  situation  occurs  since  the  trailing 
vortices  do  not  act  to  restrain  the  bound  circulation  at  this  segment, 
but  rather  become  a  destabilizing  influence  in  the  solution  for  the 
unknown  bound  circulation  strengths. 

Previously,  Queijo  has  shown  [24]  that  a  better  representation  of 
spanwise-flow  effects  on  a  sideslipping  wing  consists  of  trailing 
vortex  elements  that  remain  fixed  in  the  wing  (parallel  to  the  chord- 
wise  direction)  until  they  stream  free  at  the  trailing  edge,  a  geometry 
originally  proposed  by  Weissinger  [25].  Experience  with  the  TBL 
computational  procedure  indicates  that  modification  of  the  TBL  model 
to  incorporate  the  trailing -vortex  geometry  first  proposed  by  Weissinger 
[25]  will  have  a  twofold  benefit.  First,  the  existing  restriction  on 
advance  ratio  and  segment  widths  will  be  removed.  Second,  the 
approximation  introduced  by  dividing  the  blade  into  a  finite  number  of 
spanwise  segments  will  be  improved,  since  the  problem  associated 
with  calculating  induced  velocity  at  points  exceedingly  close  to  a  vortex 
element  will  be  eliminated. 

A  final  comment  with  respect  to  the  adequacy  of  the  adopted  wake 
model  bears  on  the  question  of  the  circulation  strengths  assigned  to  the 
rolled-up  root  and  tip  vortex.  Calculations  made  during  the  wake- 
simplification  study,  for  one  specific  rotor  and  flight  configuration, 
indicated  that  the  maximum  circulation  on  the  blade  did  not  move 
radially  to  any  significant  extent.  Further,  the  circulation  distribution 
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tended  to  be  relatively  flat  near  the  peak.  Accordingly,  it  appeared 
reasonable  to  assign  the  circulation  strength  computed  at  the  specific 
radial  station  bracketing  the  average  location  for  peak  gamma,  r  ,  to 
the  rolled-up  root  and  tip  vortices.  This  procedure  obviated  the  need 
for  program  lo^ic  to  seek  the  maximum  value  of  r  along  the  blade  span 
and  thereby  simplified  the  computational  procedure.  Calculations  made 
for  the  purpose  of  comparing  theory  with  data  obtained  in  the  NASA-Ames 
Tunnel  have  shown,  however,  that  "flight"  conditions  can  exist  wherein 
the  location  of  peak  r  moves  almost  from  root  to  tip  in  an  erratic  fashion 
during  one  rotor  revolution.  It  appears  that  the  scheme  adopted  herein 
for  assigning  circulation  strengths  to  the  rolled-up  wake  is  not  always 
a  good  approximation.  Thus,  an  error  of  uncertain  magnitude  has  been 
introduced  into  the  transient-response  computations  performed  during 
this  investigation.  Future  modeling  of  rotor-wake  systems  should 
make  due  allowance  for  the  fact  that  the  spanwise  location  of  maximum 
circulation  is  not  as  stationary  as  was  assumed  in  this  study. 
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Figure  I.  PICTORIAL  EXAMPLE  OF  THE  INITIAL  PORTION  OF  THE  WAKE  OF  A 
TWO-BLADE  ROTOR  DIVIDED  INTO  FOUR  RADIAL  SEGMENTS. 


Figure  2.  WAKE  MODELS  EXAMINED  IN  A  WAKE-SIMPLIFICATION  STUDY. 
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AERODYNAMIC  L0A0IN6,  T,  Ibs/ft 


EXPERIMENTAL  DATA  OBTAINED  WITH  HU-lA  HELICOPTER  [|8J 


O  ORIGINAL  CALCULATION:  COMPLETE  GRID  OF  TRAILING 
^  AND  SHED  VORTICJTY 


SHED  VORTICITY  IN  WAKE  TRUNCATED  AT  1/4  REVOLUTION 
BEHIND  BLADE 


SHED  VORTICITY  IS  COMPLETELY  ELIMINATED;  GRID  OF 
TRAILING  VORTICITY  RETAINED 


SHED  VORTICITY  TRUNCATED  AT  1/4  REVOLUTION 

SLADE:  TRAILING  VORTICITY  ROLLED  UP  INTO  T 
ND  TP  VORTEX 


BEHIND 
E  ROOT 


RAILING  VORTICITY  ROLLED  UP  INTO  TH 
VORTEX  l/4  REVOLUTION  BEHIND  BLADE 


SHED  VORTICITY  TRUNCATED  15°  BEHIND  BLADE:  TRAILING 
VORTICITY  "ROLLED  UP"  INTO  TIP  VORTEX  ONLV 


0  GO  ISO  ISO  240  100  MO 


AZIMUTH  ANGLE,  f  ,  digrui 


Figure  3.  INFLUENCE  OF  WAKE  REPRESENTATION  ON  COMPUTED 
AERODYNAMIC  LOADING  vs.  AZIMUTH  ANGLE. 
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AERODYNAMIC  LOADINO,  Y,  lb*/ft 


Figure  S.  SCHEMATIC  DIAGRAM  OF  TBL  MODEL. 
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Figure  6.  AXIS  SYSTEM  USED  IN  TBL  ANALYSIS. 
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Figure  7.  STEADY-STATE  AERODYNAMIC  LOADINGS  v$.  AZIMUTH  ANGLE 
H-3N  ROTOR,  /*  =  0.  18. 


PD 


r/»  =  0.25 


Figure  8.  COMPARISON  OF  INDUCED VELOCITY  DISTRIBUTIONS 
YIELDED  BY  SSBL  AND  TBL  CALCULATIONS. 


INMCI9  VELOCITY,  m,  ft/stc 
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0*1 


r/R  =  0.95 


T  -  AERODYNAMIC  LOADING,  lbt/ft 
t  -  FLAP  ANGLE,  radiant 
«  -  WAKE  TRANSPORT  VELOCITY,  ft/aac 
90-  COLLECTIVE  PITCH,  radiant 

-  MEASURED  RESPONSE 

.  CALCULATED  RESPONSE 


200  f  r/R  =  0.55 

'T  .  W  fiA 


100  Y  y  \ 


V  \  :/  v  A 
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T  .[  /A 
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AZIMUTH  ANGLE,  V  ,  dagraat 


r/R  *  0.85 
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AZIMUTH  ANGLE,  *  ,  dagraat 


Figure  10.  MEASURED  AND  CALCULATED  TRANSIENT  FLAPPING  AND  AIR  LOAD 
RESPONSES;  H-34  ROTCR,  WIND-TUNNEL  RUN  I,  fju  =  0.202. 


r\ 


MEASURED  RESPONSE 
CALCULATED  RESPONSE 


T  -  AERODYNAMIC  LOADING,  Tbs/ft 
A  -  FLAP  ANGLE,  radiant 
w  -  WAKE  TRANSPORT  VELOCITY,  ft/aac 
COLLECTIVE  PITCH,  radiana 


AZIMUTH  ANGIE,  ^dagraaa  AZIMUTH  ANGLE,  f  ,  dagraaa 

Figure  12.  MEASURED  AND  CALCULATED  TRANSIENT  FLAPPING  AND  AIR  LOAD 
RESPONSES;  H-3Y  ROTOR,  WIND-TUNNEL  RUN  3,  /*,  =  0.287. 
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360 


SPANWISE  POSITION,  r/R,  NO NO I MENS  I ORAL 

Figure  II.  INDUCED  VELOCITY  vs.  SPAN;  COMPUTED  FOR  TV©  DIFFERENT  DISTRIBUTIONS 
OF  VORTICITY  IN  THE  WAKE;  H-3N  ROTOR,  u.  =  0.18. 


AZIMUTH  AMLE,  <f  ,  dtjrni  AZIMUTH  AN  ALE,  V  ,  degree*  AZIMUTH  AN  ALE,  f  ,  degree* 

Figure  15.  TRANSIENT  RESPONSE  TO  A  RAPID  INCREASE  IN  COLLECTIVE  PITCH; 

CALCULATED  FOR  THE  H-34  ROTOR  IN  HOVERING  FLIGHT. 


0  ISO  0  ISO  0  ISO  o  ISO  O'"  ■  ’o  ISO  0  ISO  0  ISO  0  II 

AZIMUTH  ANGLE,  f  ,  digrNi  AZIMUTH  ANGLE,  V  ,  dtgrtti 


Figure  16.  TRANSIENT  RESPONSE  TO  A  RAPID  DECREASE  IN  COLLECTIVE  PITCH 
CALCULATED  FOR  THE  H- 34  ROTOR  IN  HOVERING  FLIGHT. 


TWO  BLADED  ROTOR  DIVIDED  INTO  FOUR  RADIAL  SEGMENTS 


Figure  17.  INDEXING  PROCEDURE  USED  TO  DESIGNATE  BLADE  AND  WAKE  GEOMETRY. 
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S|TJJ£  JOIN.  PROGJAM  JAJLJJ: 
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I 

© 

WRITE  HEADING  PAGE 
ON  OUTPUT  TAPE 

NTRAN8 


PROGRAM  READY  FOR 

■mar, 


SET  UP  ADDITIONAL  TABLES 
AND  PROSRM  CONTROLS: 


iiHSSSun. 

- I - 


FNR 


INITIALIZE  INDUCED  VELOCITY: 


W(J,I)  x  WIN(J) 

( (J*l  »LCN)  ,1  =  I  »NB) 


MIL  COHTHE  \ 

\amr 


/  CALL  OUTPUT  A 

^  ,t0,l  J 


Figure  16(e).  CONDENSED  FLOW  CHARTS: 


-  "BLADES". 


46 


1L/W) 
IT 


*MJFT_TA*LEJ 

^nr, 
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UfDATE_T*ILE_S 

CNR 


I 


COMPUTE.  CUMEliT 

*  I, 

I 


*nr.  *nr 


STORE  CURRENT  W 

wmTlwnT=wc“ 

I 


£0MFy  LEJJ.*0_EJ  EOM  E  TRY 
x.y.z  FOR  ALL  (J.l) 


mTLTAJLEJ 

»R.  *R 


YES 


■*l 


SOMRUTE_CyRREIIT 
*R,  ^R 
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CALL  OUTFUT 

t  YES 

0,x.Y,< 

(  <'1**)nr 

f 

(<•*!•<  )d 

<3> 


SHIFT  TABLES 

SET  CONTROLS  FOR 
ROLLED-UF  WAKE 

f.  ■ 

L - -  '1 - 

i 

SET  CONTROLS  FOR 
NONiROLLEO-UR  MAKE 


CONFUTE  CONTRIBUTION  OF 
_ _ROLLED-UP_WAjfE___ 

(  0  CHECK  FOR  TIMER  OVERFLOW 

(2  CONFUTE:  CT1^ 

ib,«,i 

NILN+la  KSNA.NR+ 1 
(S)  SUM:  W^f (rR,«r) 

( N)  IF  DESIRED:  OUTFUT  r 


I 

i 


CONFUTE  CONTRIBUTION  OF 
NONROLLED:yF_WAlE_  _ 

(i)  cm eck  "foiTtTme  OVERFLOW 
<*>  COMPUTE:  <r;M  , 

22S2MBLM 

(I)  IF:  M2,  SUN:  W^f (/*„',<?-) 

IF:  M2,  SAVE:  **FORr ITERATION 
(l)  IF  DESIREO:  OUTFUT «r 


Figurt  18(b).  CONDENSED  FLOW  CHANTS:  MAIN  PROGRAM  -  "BLADES”. 
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COMPUTE  MHTIALJAMMA 

(1)  COMPUTE:  Up.Uj  AND  SAVE 

FOR  PROCED  SUB R. 

(2)  COMPUTE:  /^(uJ.Uj.B) 


NTROL^ 


SHIFT: 
Tr  TABLE 


PREPARE  MATRIX  FOR  SIMSOL 


(1)  INPUT:  «*FOR  1=2 

(2)  COMPUTE:  «rT=f (BJ.  <r,i ,J , ib.C) 


SOLVE  FOR  CURRENT  BAMMA 


( 1)  CHECK  FOR  TIMER  OVERFLOW 

(2)  INITIALIZE  rNLaf(aT’,J’UT> 

(3)  PREPARE  SIMSOL:  rg-^L 

(»)  CALL  SIMSOL  FOR  f] ’M 

(6)  SUM:  W^ft'M’^o) 

(6)  COMPUTE  CURRENT  rNL  =f(arT,BJ,uj,wA) 

(•)  IF  ANY  CTj  >0,  PREPARE  FOR 
PERCENTAGE  TEST  ON  C 

(b)  IF  PERCENT  TEST  FAILS 
1  RETURN  TO  (3) 

(c)  IF  ALL  TESTS  PASSEO,  fj.  IS 
FINAL  SOLUTION  AND  W=f(WT,WA) 


SHIFT: 
TNr  TABLE 
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*,(•). 1=1, NB 


rwr 

TIMER 

^OVERFLOW, 
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JO  STEADY- STATE 
TEST  TABLE 


^steaoyN 

STATE 
L ACHIEVED; 


SAVE  ON  TAPE  I 
lREQU I  RED  INFORMATION 
POf  PLOTTING  PROGRAM 


TO.W*,, 


Figuro  18(c).  CONDENSED  FLOW  CHARTS:  MAIN  PROGRAM  -  "BLADES". 
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ENTRY 


FRESET  PATH  CONTROLS 

AND  TABLE  INDICES 

irfr 

©A 

> 

ZERO  OUT  SUNS  AND  SET 
RUNGE-KUTTA  CONTROLS 

'.u.  =  0 
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=  0 
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*2  =  ^T 
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SET  TABLE  INDEX  FOR 
VARIABLES  CC,  UT,  /*T 
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CONFUTE:  (Ojj,  j= I.LCN) 
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,(V 


i 

AND  6B.  MIEN  AFFLI  CABLE 
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\  "° 

OF  RUDE  i  FOR  USE 

IN  CONFUTATION  OF 

/ 
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l 


CALL  FROCED  SUBROUTINE: 

( 1)  CONFUTES  QT.Uj,  f*T  MIEN  NEEDED 
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TAFE  MIEN  REQUIRED 


CONFUTE: 

•*  =  f  *1 ) 

1 — 


CONFUTE: 
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Fig-rt  19(a).  CONDENSED  FLOW  CHARTS:  SUBROUTINE  -  "COMEQU". 
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STOR(.  rot  OUTPUT 
fi\  =  if 


COMPUTE: 

(1)  Ahp.Au^  *  f(*0**» 

(2)  lhF  =  f(f (»♦*), !(•♦*)) 


AND  40  WHEN  APPUCAILE 


Figurt  19(b).  CONDENSED  FLOW 


Figure  20.  CONDENSED  FLOW  CHART:  SUBROUTINE  -  "PROCED". 
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BIT  IT 
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COMPUTE  WAKE  CONTRIBUTION 
WC  = 


NO 


RESET  fi  CHECK  PATH,  1SSA  «  4-1 
SET  AW  BYPASS  SWITCH,  ZSWI  =  0 


Figure  22.  CONDENSED  FLOW  CHART:  SUBROUTINE  -  "COMWB". 
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APPENDIX 


THE  MATHEMATICAL  MODEL  AND  COMPUTATIONAL 
PROCEDURE  USED  TO  CALCULATE  TRANSIENT 

BLADE  LOADS 


INDEXING  NOMENCLATURE 

A  concise  format  for  the  equations  forming  the  mathematical  analog  of 
a  multibladed  rotor  generating  a  vortical  wake  requires  the  use  of  six 
indices.  These  indices  are  used  both  as  subscripts  and  superscripts 
with  the  variables  in  the  equations.  Their  purpose  is  to  designate  the 
time  and  spatial  dependence  of  each  variable  under  consideration. 

The  indices  are  m ,  L  ,  j  ,  ib  ,  m  ,  and  J.  They  designate,  respectively 
(1)  an  instant  of  time,  (2)  a  blade  of  the  rotor,  (3)  a  spanwise  location 
on  a  blade  (corresponding  to  the  midpoint  of  a  designated  segment),  (4) 
the  wake  trailed  by  blade  i  ,  (5)  an  azimuthal  position  in  the  wake  and 
(6)  a  spanwise  position  in  the  wake  (corresponding  to  the  endpoints  of 
the  segments  into  which  the  blades  have  oeen  divided).  Figure  17 
pictures  the  manner  in  which  these  indices  are  used  to  designate  blade 
and  wake  geometry,  for  example.  Although  the  time  index,  *?,  does 
not  appear  on  the  Figure,  it  should  be  understood  that  a  value  of  mis 
associated  with  the  azimuthal  position  of  a  blade  designated  as  the 
reference  blade  (usually  blade  ('  :  1)  since  time  and  azimuth  are 
equivalent  variables  for  a  rotor  rotating  at  a  fixed  r.  p.  m. 

INPUT  DATA 

The  input  data  required  to  compute  transient  air  loadings  fall  into  the 
following  categories: 

(1)  Data  specifying  the  geometry  of  the  rotor  and  its  wake, 
including  the  specification  of  the  number  of  calculation 
points. 
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(2)  Data  specifying  the  equilibrium  flight  condition,  including 
the  trim  values  of  collective  and  cyclic  pitch. 

(3)  Data  specifying  the  mass  properties  of  the  rotor. 

(4)  Starting  values  for  flap  angle,  flapping  angular  velocity, 
and  induced  velocities. 

(5)  Data  specifying  the  variation  in  collective  pitch  and 
in  the  mean,  induced  downwash  with  time. 

(6)  Data  necessary  to  control  the  operational  logic  built  into 
the  digital -computing  program. 


WAKE-GEOMETRY  COMPUTATIONS 


The  coordinates  of  the  midpoints  of  the  blade  segments  (at  the  three- 
quarter  chord)  in  the  adopted  i-y-f  axis  system  (Figure  6)  are 
designated  as  [*m;j  l  </m  j  j  i  fm  ^  ] •  Similarly,  the  coordinates 
of  the  endpoints  of  the  elements  of  trailing  vorticity,  in  this  same 
axis  system,  are  designated,  as  [  ]• 

The  azimuthal  position  in  the  wake  at  which  the  wake  is  assumed  to  roll 
up  into  two  vortices  (i.  e. ,  one  at  the  blade  tip  and  one  at  blade  root)  is 
designated  as  m  -  MBLM.  With  this  notation,  it  becomes  convenient  to 
write  a  separate  set  of  equations  defining  the  geometry  of  the  nonrolled 
wake  ( HittfSltf)  and  the  rolled -up  wake  (W2  MBU4), 


With  the  aid  of  the  definitions  of  and  Fg  (given  in  Figure  17)  and  of 
the  nomenclature  called  out  in  Figure  6,  and  on  further  defining 

W 'm  -  mean  induced  downwash  governing  the  vertical  transport 
of  the  vorticity  in  the  wake  (W  is  a  function  of  time) 


?(J"t,2)  =  radial  distance  to  spanwise  location  of  rolled-up  root 
and  tip  vortex, 

we  obtain  the  following  set  of  equations  for  [  t,  y,  f  ]  and  [£,  q ,  C  ]: 


X,  u,  ]  of  Blade-Segment  Midpoints  (at  3/4  chord) 


"  rjco*  *  bjsin  .  -  (rj  -€ )  <*g/8m  ^ 


im.i.j  *  'i  *" +m.i  '  *•*  *m.i 


*  n<*.  c°*  * *.<*, «» 


[flai 


of  Vortex-Segment  Endpoints;  Nonrolled  Wake 


iii  a  MB  LM  ;  m  <  m 


m  *  ’  ■  im.it.  m.i  -(*-•>  «.  A.,.* 


m  *  *  '  im.it, m.i  '  im-t.it.  m-t.i 


CL(NA) 


m-1:  7 •  rJ'sin 


m  m/  ;  tm.ik.rn, A  *  cos 


2  TT 

m>f  '  tm.Lk.m.A  *  *  £1(NA)  "m-m-t 
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f  ]  of  Vortex-Segment  Endpoint*;  Rolled-Up  Wake 


M3L  M  S  m  £  NA  (NR)  +  1  ;  m  <  m 


**mN3LM’  'tm.ik'MH.Hti  *  O  c°$iPm-MSLM+l,ib  “  ^  r4  ~g)  asAm-M»U**1.ih 

+(MeL"-'&)V 


m>M3LM .*  j  +  a(#A)  V 


9Ft  »  MBLM:  *44*.  J  ^  ~€)Am-  MBiH+t,  *»  +  “j  as  COS  +m.M0L*+,,U 


nw*;  *»- 


,  JL  *  Cn</,  4^  S-f  </ 

where: 


2tt  _ 


f  # 


^m,ih  "  »  ^*,4*  *  Am.i 


-*•**+ 


f,.'4  *  *m-*9LM+t ,i 


Am-*  +  t,  ik  *  fim-m*  f,  4 


Am-M$LM+ 1,  tS  "  Am  -  *44  M+1,i 
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COMPUTATION  OF  a-  MATRIX 


The  biot-Savart  law  yields  the  relationship 

w  •  a-y 

for  the  velocity  induced  in  the  }  direction  by  an  arbitrarily  oriented 
straight  vortex  filament  of  constant  vortex  strength,  y.  The  coefficient 
cr  is  a  function  only  of  the  coordinates  of  the  point  where  the  velocity 
is  being  computed  and  the  coordinates  of  the  vortex  filament  endpoints, 
i.  e. , 

A  complete  set  of  <r  coefficients  represents  the  influence  of  all  of  the 
elements  of  vorticity  in  the  wake  (of  unit  strength)  at  all  of  the  load 
calculation  points  on  each  blade. 

It  is  convenient  to  indicate  the  coefficient,  <r  ,  applicable  to  a  specific 
element  of  trailing  vorticity  and  to  a  specific  segment  calculation  point 
as  .  The  superscripts,  m,  i  ,  and  j  ,  indicate  the  specific 

location,  as  given  by  its  [sm  iJ  i  ]  coordinate,  at 

which  the  induced  velocity  is  being  determined.  The  subscripts,  ik  ,  m  , 
and  J  indicate  the  element  of  vorticity  in  the  wake,  as  given  by  the 
f  ,m,J  ;  7 ;  Cm.tk.TK'il  coordinates  at  both  the  fore  and 
aft  ends  of  the  element.  The  notation  is  such  that  the  aft  end  of  the 
vortex  element  is  located  at  the  azimuthal  position,  m ,  with  the  forward 
end  at  the  azimuthal  position,  m  -  1.  Accordingly,  no  value  can  exist 
£or  when*s  '• 

For  me 2,  application  of  the  Biot-Savart  law  yield 
a.m,4  *  *-4-4  *  m^mk 

where 
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m+  * 

*  Jm, . ,  j  “  Gm,  4 


mk  *  y**.*.j  *  *1**,  ik.  ifi-  t,  J 
nb  *  ~ 

rj  -•  ^  *  V 

*  j  I  .  S  g 

•  U*  +  "k  + 


The  above  formulation  for  o-  is  unnecessarily  lengthly  for  vortex 
elements  located  at  a  large  distance  from  a  blade  segment.  Thus  the 
above  equation  is  used  for  rj<  ,  where  rmf  is  an  input  quantity  to 
the  computer  program.  For 

■^mmb  m«» 


r~,  2  r*  5  r~«  • 


firrJ 


*£>  rmt  : 

INDUCED-VELOCITY  COMPUTATION 


'  * 


For  purposes  of  implementing  the  TBL  computation,  the  contributions 
of  all  vortex  elements  to  the  total  velocity  induced  at  a  blade-segment 
location  are  summed  over  the  entire  wake  with  the  exception  of  the 
elements  located  immediately  adjacent  to  the  blades.  It  should  be 
noted  that  the  vortex  strengths  of  the  latter  elements  are  a  function  of 
the  bound -vorticity  distribution  on  the  blade  and,  thus  are  unknown 
quantities  prior  to  the  determination  of  these  bound  vorticities. 

Rather  than  require  the  program  to  search  for  the  maximum  value  of 
circulation  on  the  blade  for  eventual  assignment  of  these  strengths  to 
the  rolled-up  root  and  tip  vortices,  we  arbitrarily  specify  that  the 
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circulations  associated  with  a  given  spanwise  segment  be  assigned  to 
the  rolled-up  wake.  Past  experience  is  employed  for  guidance  in 
selecting  that  spanwise  segment  which,  on  the  average,  sees  the 
maximum  circulation  along  the  span.  We  designate  the  incomplete 
summation  for  the  induced  velocity  as  wr  and,  in  accordance  with  the 
above  remarks,  express  wr  as: 

mb  Cmblh  ItN+t  , 

S  *at 


where: 


MA(MR)r1  .  .  n  1 

*  H  rik.  m.m  +  i,  JSO  -  Cri?,M'U J  | 


m  •  MiLM+t 


m  2  2 


and 

m  -tB  +2 ,  £  8 

with 

/**'•  m-fR  +  2 ,  i  ■  ib ,  i=U> 

with 

m'  •  m  -  m  *-2y  i  •  ib  ,  J  »-/  -  / 

ib.m-  tR  +  Z,  J 30 

with 

m'w  m  -m  *-2,  i  •  ib  ,  J  *  JSG 

N6 

=  number  of  blades 

LCN  =  nvunber  of  blade  load  points 
NR  -  number  of  revolutions  of  wake 

J5<r  =  blade  station  associated  with  maximum  circulation. 


ROTOR  KINEMATICS 

For  a  given  blade  segment,  J  ,  located  at  a  distance,  r-  ,  from  the 
center  of  rotation,  the  aerodynamic  angle  of  attack,  measured  in  a 
plane  perpendicular  to  the  blade  span,  can  be  expressed  as 

/T  9  0  ^ 

where 

0m  ■  j  is  the  geometric  pitch  angle  of  the  blade  section 
measured  with  respect  to  the  shaft  plane 

and 
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'Up  .  . 

•  •  *  tan’  "**'«* 

^  **  «/  U- 

**  •  -  1  • 

A 

The  "perpendicular"  and  "tangential"  components  of  the  airstream 
relative  to  the  blade,  at  the  three-quarter  chord,  are: 

UPm  i  J  c  ~Va*  *  "m.i.j  '“m.i.j  **’[V+  "m.L.j  <*»  +  “m.i.j]  fim.i  COS  *m,i 

mV,mtijArn.i  *  bi  {**m’n  [?tc*n  ’  °*  Co*  j 

urm,ij  •  [v*  «,  «■  <L.i  ‘  cos  *m.L  *  fl  r, 

where 

unt  i  Ji  i  J  are  ^e  }  components,  respectively, 

of  the  velocity  induced  by  the  wake  at  time,  m,  at 
blade,  i  ,  and  station,  j . 

00  is  the  rate  of  change  of  collective  pitch  (a  function 


of  time) 


and 


0tc  and  9ta  are  the  cosine  and  sine  components  of  cyclic  pitch, 
respectively. 


Otg  and  fi  have  been  considered  to  be  small  angles. 


It  should  be  noted,  that  whereas  some  of  the  terms  in  the  equation  for 
Up  may  be  significant  at  very  low  advance  ratios,  at  moderate  flight 
speeds  we  have 

a,  v*  <<  V 
w  <x$  «  V 
ujttvy8  <<  V*9,w 
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The  equations  for  Up  and  Uj  thus  simplify  to  the  following  set: 

’  'V“*  *  i 

rbi  -  6i.  «"  *<-,.]) 

Of  .  .  »  Vsut  4>  •  *  fir- 

As  was  done  in  Reference  12,  the  velocity  component,  6^  ,  is 
considered  to  consist  of  (1)  the  velocity  induced  by  the  wake,  i.e. ,  w , 
and  (2)  the  velocity  caused  by  the  "effective  plunging"  velocity  of  the 
blade.  Thus  we  have  that 


Up  .  .  ■  tv  .  .  -  h  •  . 

4  «  J  i  |  6  0  J 


where 


-  tj  { e,m  -  n  [«9/e  3*>i  ^  cos  7\] 

Finally,  we  note  that  the  geometric  pitch  angle,  j  ,  of  a  given 
blade  segment  is  given  by  the  following  kinematic  relationship: 


where 


j  •  e.m+—  fi  *  9k  «•*  *m,i  *  *m.i 


0o  is  the  collective  pitch  setting  at  the  blade  root  (or 

m 

at  any  specified  spanwise  location)  and 


is  the  twist  (assumed  linear)  of  the  blade  (rad/ft) 
SOLUTION  PROCEDURE  FOR  UNKNOWN  CIRCULATION  STRENGTHS 


As  was  discussed  in  the  body  of  the  report,  the  determination  of  the 
unknown  bound  vorticities,  rm  -  •  ,  is  facilitated  by  linearizing  the 
rotor  kinematics,  namely  by  assuming  that 

,  Uhm  i  ;  UP  •  • 

•  «  tm*  ' 

(Jm  •  .  (/*  •  • 
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On  so  doing,  we  have  that 


where 


U 


0 


U, 


I  •  .  +  Um  .  1 


a*  s 


dCt 

<£ac 


(considered  as  a  constant  with  no  m,  *  ,  or  J 
dependency). 


Since  Up  .  .  is  a  function  of  **'  ;  •  which,  in  turn,  is  a  function  of 
all  the  rm  (i.e. ,  for  all  i  and  j  ),  the  above  equation  contains,  in  an 
implicit  fashion,  all  of  the  unknown  circulation  strengths  existing  on  the 
blades  at  time,  m  .  A  set  of  equations  relating  each  rm  ;  ;  to  all  of  the 
rm^  j  's  existing  ^t  time  m,  must  thereby  be  solved  simultaneously  to 
determine  the  unknown  circulation  strengths. 


In  order  to  facilitate  this  simultaneous  solution,  we  define  several  of 
the  existing  variables  or  quantities  in  terms  of  the  accumulated  contri¬ 
bution  of  the  wake  with  the  exception  of  the  trailing  elements  adjacent 
to  the  blades.  Thus,  we  employ  the  previously  defined  incomplete 
summation  for  the  induced  velocity  to  obtain  the  following  expression  for 
the  bound  circulations  due  to  the  following  known  quantities: 


Jm.i.j 


bj  0/ 


II  Ur  ■  •  I  Or  '  +U’.  • 


where 


m,i,S  • 


The  equation  for  rm  •  j  can  now  be  written  as 

#*  ICM 

*"’*J  "  a,  £  £  r j 

If  we  further  define 


iS,4  *  K+  | rik.t.4+i\ 
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where 


4  i.j 


and 


(T-V 


i  t  -  I  /M  *  tj  ^9t$4  | 

'->J»\  J.J  -  *t*mg.£94  +  fmj  +  fJ 


the  iuj  set  of  equations  for  /*_  •  •  can  be  written  in  matrix  notation 


as 


The  desired  solution  can  be  expressed  as 


This  solution  procedure  (*.s  outlined  above)  holds  for  a  linear 
relationship  between  lift  and  angle  of  attack  wherein  no  accounting  is 
made  for  the  influence  or  effects  of  aerodynamic  stall.  It  is 
possible,  however,  to  introduce  the  effect  of  stall  in  an  approximate 
manner.  In  this  study,  it  was  assumed  that  «  =  5.73  and  that  the 
blade  segments  stall  at  fourteen  degrees  with  the  lift  and  circulation 
remaining  equal  to  the  value  existing  at  ot-  14  degrees  at  any  angle 
of  attack  in  excess  of  fourteen  degrees. 


Bf'und  circulations  were  limited  to  a  maximum  value  by  computing 
the  angle  of  attack  of  each  blade  segment  and  introducing  a  limiting 
function  on  ffm  .  .  .  The  blade  segment  angles  of  attack  are: 


Ui  #A  .  . 

oc  •  .  ■  0  -  •  + 

m.t.j  JJ1  — 

where  the  velocities  induced  by  trailing  elements  adjacent  to  the 
blades  are  designated  as  wA  .  .  ,  and  are  given  by 


Note  that  rm  ^  j  are  unknown  quantities  until  the  equation  set  has 
been  solved  at  least  once.  Accordingly,  the  are  first 

evaluated  with  .  .  assumed  to  be  zero.  If  any  of  the  initially 
determined  *m,ij  exceed  the  specified  stall  angle,  we  solve  an 
equation  set  of  the  form 


where 


I |  a(l  I  "  asUll| 

subject  to  the  requirement  that 


Only  one  solution  of  the  simultaneous  equation  set  is  required  if  all 
the  .  .  are  equal  to  zero.  If  this  quantity  exists,  the  equation 

set  is  solved  and  the  .  .  are  redetermined.  The  solution 

procedure  is  repeated  until  the  changes  computed  for  Cm  •  •  are  less 
than  a  specified  percentage. 


BLADE -FLAPPING  DYNAMICS 


The  aerodynamic  loading  is  designated  as  Tm  •  •  and  under  the 
assumptions  made  in  the  study  can  be  expressed  as 

7*  .  •  *  I®  Ot  •  •  r* 

The  moment  about  the  flapping  hinge  caused  by  the  aerodynamic 
loading,  when  expressed  as  a  continuous  function,  is 

*Vi  ’  T(r)dr 

and  for  the  purposes  of  this  study  was  numerically  evaluated  by  means 
of  a  trapezoidal  integration  procedure.  The  resulting  aerodynamic 
moments  are  used  in  the  blade. flapping  equation: 
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to  obtain  the  angular. flapping  accelerations  at  each  instant  of  time. 

It  should  be  noted  that  the  Runge-Kutta  integration  procedure,  used  to 

• «  • 

integrate  4  to  yield  and  ,  utilized  predictions  over 

the  half  and  full  interval  with  WTmt  •  ;  assumed  to  be  constant  over 
the  time  (azimuth)  interval.  No  attempt  was  made  to  include  the  wake, 
geometry  and  induced. velocity  computations  into  the  Runge.Kutta 
integration  since  this  procedure  would  have  been  prohibitive  in  terms 
of  computing  time  and  cost.  Since  the  coupling  between  Wr  and  blade 
flapping  during  a  single  time  interval  is  very  small,  the  errors  intro- 
duced  by  removing  the  Wj  calculation  from  the  Runge-Kutta  procedure 
are  believed  to  be  negligible. 

DIGITAL-COMPUTER  PROGRAM 

The  above  outlined  mathematical  model  was  programmed  by  means  of 
the  IBM  FORTRAN  IV  language  for  the  7044  computer.  The  limited 
storage  capacity  of  the  7044  required  that  a  small  portion  of  the 
coefficients  be  computed  and  then  used  immediately  to  compute  their 
contribution  to  Wr  .  Since  the  e* coefficients  for  the  elements  adjacent 
to  the  blades  must  be  saved  for  use  in  the  simultaneous  solution 
procedure,  lack  of  sufficient  storage  capacity  required  that  these 
coefficients  be  read  out  on  tape.  Although  this  requirement  has  been 
a  significant  factor  in  determining  the  computing  time  per  azimuth 
interval,  the  major  time-consuming  portion  of  the  program  is  the 
wake  geometry  and  ^-calculation  phase.  For  the  case  of  a  four- 
bladed  rotor  with  three  revolutions  of  wake,  the  program  in  its 
present  form,  requires  approximately  0,64  minutes -per -azimuth 
interval  when  no  iterations  are  required.  Since  100  time  intervals  is 
a  typical  example  of. the  number  of  intervals  required  to  compute  either 
a  steady- state  or  transient  response,  a  typical  computation  requires 
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approximately  one  hour 


Condensed  flow  charts  for  the  main  program  and  the  subroutines  are 
presented  in  Figures  18  through  22.  These  flow  charts  can  be  compared 
with  Figure  5  to  show  the  manner  in  which  the  developed  mathematical 
analog  was  augmented  with  additional  logic  to  facilitate  numerical 
computations. 
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